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rile  loading  data,  obtained  from  gas  gun  experiments,  was  fitted  by  the  function 


Cr  (kbar)  = 13.65 


m- 


4.358  ^ 


• 9227(f)  . 0.0684(f)" 


where  a is  the  compressive  stress  in  the  direction  of  wave  propagation  and  AK/R  is 
the  fractional  change  in  electrical  resistance.  The  standard  deviation  associated 
with  the  above  expression  is  0.638  kbar.  Better  prediction  at  lower  stress  levels 
can  be  obtained  by  using  a different  equation: 


0-  (kbar)  = 1.082 


-[20.80(f) 


+ 9.168 


The  standard  deviation  associated  with  this  equation  in  the  range  0-10  kbar  is  0.252 
_kbar.  The  ne^d  for  the  two  equations  arises  from  the  shape  of  the  loading  curve  in 
the  1 to  2 kb  r region— the  change  in  slope  which  occurs  is  not  adequately  described 
by  the  fourth  order  polynomial  that  fits  well  over  the  entire  range. 

> 

Between  0 and  about  10  kbar  the  unloading  paths  are  a family  of  straight  lines 
with  both  slope  and  intercept  being  functions  of  peak  stress.  The  release  data 
from  stress  levels  between  about  10  and  33.40  kbar  have  been  fitted  by  second  and 
third  degree  polynomials.  Stress  cycling  affects  the  response,  but  insufficient 
data  were  obtained  to  quantitati vely  describe  reloading  behavior. 

The  maximum  stress  at  which  ytterbium  can  be  used  safely  is  approximately  30 
kbar  because  of  the  occurrence  of  a polymorphic  phase  transition  between  30  and  35 
kbar  that  is  essentially  complete  in  a few  microseconds. 

An  attempt  has  been  made  to  assess  the  effects  of  the  tangential  components  of 
a divergent  flow  field  on  the  response  of  ytterbium  transducers  through  experiments 
designed  to  investigate  the  tensile  strain  dependence  of  electrical  resistance. 
Difficulty  in  ascertaining  how  much  tangential  strain  is  actually  coupled  into  a 
gage  in  a divergent  flow  field  precludes  application  of  the  data  at  this  time. 

It  was  ascertained  that  the  piezoresistance  coefficient  of  annealed  foil  was 
much  greater  than  the  coefficient  of  the  as-rolled  foil  normally  used  for  gages,  but 
that  the  greater  part  oi  the  resistance  change  during  loading  is  irreversible.  In 
addition,  it  has  been  noted  that  at  stresses  of  10  kbar  and  below  there  is  a large 
difference  betweeen  the  hydrostatic  piezoresistance  coefficient  and  the  coefficient 
measured  by  gas  gun  experiments.  The  difference  is  greatest  at  about  1,5  kbar. 
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1.  INTRODUCTION 


Prediction  of  the  iui  situ  response  of  rocks  and  soils  to  dynamic 
loading  requires  accurate  knowledge  of  the  constitutive  relations  among 
stress,  particle  velocity,  and  density  during  the  passage  of  a stress  wave. 
Good  ±_n  sit u measurements  would  be  the  best  source  of  data  from  which 
constitutive  relations  could  be  obtained  because  these  measurements 
include  the  effects  of  cracks,  voids,  and  residual  stresses.  Such 
measurements  of  stress  in  rocks  and  soils  have  been  made  with  transducers 
having  ytterbium  as  the  stress-sensitive  element.  Keough  (1970) 
first  suggested  ytterbium  for  use  in  this  application,  primarily 
because  it  exhibits  a very  large  positive  change  in  electrical 
resistance  with  compressive  stress,  and  the  ratio  of  this  signal  to  any 
divergent  flow-induced  strain  signal  was  expected  to  be  small  over  the 
stress  range  of  interest.  Although  high  quality  resistance-time  records 
have  been  obtained  in  the  field  during  a number  of  nuclear  and  high  explosive 
detonations  (for  example,  see  Smith,  1972;  Grady  et  al.,1973),  conversion 
to  stress-time  histories  has  been  hampered  by  a lack  of  data  on  the  relation- 
ship between  stress  and  resistance  for  ytterbium. 

The  program  described  in  this  report  was  designed  primarily  to 
increase  the  utility  and  credibility  of  ytterbium  gages  for  in  situ  rock 
and  soil  stress  measurements.  However,  the  calibration  data  and  other 
results  are  applicable  to  ytterbium  gages  used  in  other  situations.  The 
specific  objectives  were 

(1)  To  generate  loading  and  release  calibration  data 
immediately  applicable  to  the  analysis  of  gage 
records  now  being  obtained  in  the  field  and 
laboratory  with  ytterbium  stress  transducers. 

(2)  To  improve  the  accuracy  and  reliability  of  ytterbium 
stress  transducers  through  investigation  of  some  of  the 
existing  questions  about  the  piezoresistance  of  ytterbium. 
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Calibrating  dynamic  stress  gages  oi  the  type  used  in  this  program 
requires  the  following  experimental  methodologv.  The  material  to  be 
calibrated — yt terbium  in  this  instance — is  fabricated  into  some  con- 
venient form  and  sandwiched  between  two  pieces  of  a material  whose 
stress-particle  velocity  Hugoniot  has  been  ascertained  through  the  use 
of  some  other  technique,  such  as  laser  interlcromeLi y . When  we  impact 
this  target  with  a projectile  of  either  the  same  material  or  some  other 
material  with  well-de  nod  shock  properties  and  measure  the  projectile 
velocity,  we  can  spec  i>  the  peak  stress  readied  in  the  target,  Assuming 
Lite  specimen  or  gage  material  comes  to  stress  equilibrium  with  the  target, 
we  have  also  specified  the  stress  level  in  the  specimen.  Measurement  of 
the  fractional  change  in  resistance  at  peak  stress  gives  one  point  on  a 
loading  calibration  curve  at  the  stress  level  calculated  from  the  pro- 
jectile velocity  and  the  llugoniots  of  the  target  am1  projectile,  The 

locus  of  these  points  on  a s tress-, * K/H  plot  is  the  loading  calibration 

o 

curve.  We  normally  assume  that  the  gage  follows  this  curve  during  the 
loading  part  of  the  stress  cycle. 

This  program  differed  somewhat  from  other  calibration  programs  in 
that  great  emphasis  was  placed  on  obtaining  unloading  calibration  data 
so  that  complete  stress-time  histories  could  be  produced  from  gage 
records.  Experiments  designed  to  obtain  both  loading  and  unloading 
calibration  data  are  no  more  difficult  to  perform  than  experiments  that 
produce  only  loading  data,  but  the  selection  of  projectile-targe t systems 
is  restricted  by  the  need  for  detailed  knowledge  of  the  unloading  be- 
havior of  target  and  projectile  and  the  proper  impedance  mismatch. 

In  addition  to  and  concurrent  with  the  basic  calibration  phase  oi 
the  program  (which  was  of  first  priority),  we  attempted  to  gain  more 
understanding  of  the  material  parameters  that  affect  gage  behavior  by 
examining  the  gage  material  and  trying  to  correlate  differences  in,  for 
example,  metallurgical  history  with  piezorcsis tan L response. 
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I' inally,  we  approached  the  problem  of  divergent  flow  effects  on  gage 
response  by  investigating  the  effeets  of  noncompressive  strain  on  the 
electrical  resistance  of  ytterbium. 

This  program  was  not  designed  to  produce  a handbook  for  ytterbium 
gage  users;  proper  application  of  the  data  presented  in  this  report  will 
require  some  effort  and  experience.  The  production  of  a handbook  might 
be  a reasonable  activity  to  be  undertaken  in  the  near  future. 
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SUMMARY  OF  RESULTS 


The  Loll  owing  is  a summary  of  the  results  obtained  during  this  program. 
These  results  are  presented  in  more  detail  and  discussed  in  Section  1.6  and  1.7. 

2 . 1 Loading  Calibration 

* 

The  ratio  of  change  in  resistance  to  initial  resistance  (AR/R  ) 

o 

hereafter  referred  to  as  resistance  change,  of  nominal  0.005-cm  (2-mil) 
ytterbium  ‘oil  was  measured  over  the  stress  range  0.5  to  33.1  kbar,  and 
the  data  can  be  represented  in  the  range  from  0 to  10  kbar  by  the  expression 

( r AH.-)  £R 

(T  = 1.082  1 - exp  r-20.80  (— ) + 9.168  (— ) (1) 

( o ) o 

The  standard  deviation  around  the  mean  of  Lhe  curve  represented  by 

Equation  (l)  is  0.262  kbar.  A graphical  representation  is  contained  in 

Figure  33  in  Section  1.6.  There  were  33  data  points. 

We  chose  to  represent  the  entire  population  of  loading  data  (0  to 
33.10  kbar)  by  the  function 

2 3 1 

TR  Alt  AIL  ar 

Cr  = 13.65  (— ) - 1.358  (-^-)  + 0.9227  ( — ) - 0.0681  (— ) (2) 

o o o o 

The  standard  deviation  around  the  mean  ol  Equation  (2)  is  ±0.638  kbar. 

There  were  53  data  points.  The  need  lor  the  two  equations  arises  from 
the  shape  of  the  loading  curve  in  the  1 to  2 kbar  region — the  change  in 
slope  •Itich  occurs  is  not  adequately  described  by  the  fourth  order 
polynomial  that  fits  well  over  the  entire  range.  See  Section  1.7  and 
Figure  23  in  Section  1.6. 

2 . 2 Unloading  Calibration 

A number  of  unloading  (release)  paths  were  experimentally  determined 
in  the  stress  range  between  about  0.75  kbar  and  33  kbar,  Up  to 

"It  instantaneous  - initial 

K It 

o initial 
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about  1 kbar,  the  loading  path  and  the  u llcading  paths  coincide.  Between 
0 kbar  and  10  kbar,  the  unloading  paths  can  be  represented  by  a family 
of  straight  lines  of  the  form 


a = A + B (&R/R  ) 
u o u 


(3) 


where  both  A and  B are  functions  of  the  peak  stress  and  a and  (£R/R  ) 

u o u 

are  the  stress  and  resistance  change,  respectively,  along  the  unloading 
path.  Between  about  10  kbar  and  33.4  kbar,  the  measured  unloading  paths 
can  be  represented  by  curves  of  the  form 


E,  ,AR  n 

A„  ‘r’u 


i = 2 or  3 


(4) 


n=0 


o 


where  the  are  all  functions  of  the  peak  stress,  and  and  (£R/R  ) 

u o u 

are  the  stress  and  corresponding  resistance  change  along  the  unloading 
path.  See  Section  4.7.3. 

2 . 3 Reloading 

The  dynamic  piezoresistance  function  of  ytterbium  is  not  strictly 
constant  over  successive  loading  and  reloading  cycles,  but  we  do  not 
believe  the  errors  in  stress  measurement  attributable  to  this  nonconstancy 
to  be  appreciable.  See  Section  4.7.5. 

2 .4  Divergent  Flow  Effects 

A series  of  experiments  was  performed  to  determine  the  effect  of 
tensile  strain  (in  the  plane  of  the  foil)  on  the  electrical  resistance 
of  ytterbium  foil,  thereby  partially  simulating  a divergent  flow  field. 

The  resistance  shows  an  initial  decrease,  reaches  a minimum  of 
R R^  -0.02  at  about  1 . 5CL  strain  and  then  increases  linearly  with 
the  application  of  more  strain.  See  Section  4.7.8. 

2 • 5 Characterization  of  Ytterbium 

Representative  values  of  resistivity  and  hardness  were  measured, 
resistivity  was  found  as  a function  of  temperature,  and  X-ray  diffraction 
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data  were  obtained  on  the  various  types  of  ytterbium  toil  used  so  that 
some  of  the  basic  parameters  affecting  the  piezoresistance  of  ytterbium 
could  be  known.  See  Sections  1.3.1,  -1.3.2,  1.3.3,  and  -1.3.6. 

2.6  Mechanical  Properties  of  Ytterbium 

Data  were  obtained  on  the  mechanical  properties  of  ytterbium  foil 
in  tension  and  of  ytterbium  ingot  material  in  compression,  and  the 
data  were  related  to  the  observed  piezoresis tance  effects.  See  Sections 
1.3,  l and  1.3.5, 

2.7  Comparison  with  Hydrostatic  Piezoresistance  Experiments 

The  difference  between  the  piez.oresistance  under  hydrostatic 
loading  and  one-dimensional  or  two-dimensional  strain  (in  the  elastic 
region)  is  explained  by  determining  the  coefficients  in  the  piez.oresistance 
tensor.  The  analysis  was  extended  into  the  plastic  region.  See  Appendix  A. 

2 . 8 Determination  of  Upper  Stress  Limit  for  Ytterbium  Gage  Use 

The  upper  limit  tor  use  of  ytterbium  was  determined  by  shocking 
foil  gages  to  a constant  stress  above  the  pressure  required  for  the 
phase  change  under  hydrostatic  loading  and  the  change  was  found  to  occur 
at  as  low  as  about  32  kbar  under  dynamic  loading.  See  Section  4.7.6. 

2.9  Effects  of  Metallurgical  History  on  Piez.oresistance 

The  electrical  response  of  annealed  and  unannealed  (as-rolled)  foil 
material  was  compared  in  dynamic  compression  and  in  static  tension,  and 
the  state  of  hardening  was  identified  as  a major  variable.  See  Sections 
4.7.4  and  4.7.5. 

2.10  Time  Effects 

Static  experiments  were  performed  and  the  results  compared  with 
dynamic  experiments;  resistance  changes  were  measured  at  constant  loads 
over  longer  times.  No  time  effects  are  expected  for  times  in  the 
millisecond  range,  but  effects  were  observed  during  static  loading 
rates.  See  Section  4.7.7. 
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3.  CONCLUSIONS  AND  RECOMMENDATIONS 


3.1  Conclusions 

3.1.1  Calibration  Data  — In  tho  range  from  0 to  10  kbar,  the 
loading  calibration  curve  lor  as-rolled  ytterbium  foil  is  described 
adequately  by  a single  function,  Equation  (1).  The  unloading  paths 
lrom  peak  stress  up  to  about  10  kbar  can  be  said  to  form  a family  of 
straight  lines  where  the  slope  and  intercept  are  functions  of  peak  stress. 
The  loading  calibration  path  shows  curvature  above  about  10  kbar  up  to  the 
highest  stress  considered,  and  the  release  paths  become  nonlinear  as 
well.  After  attempting  to  fit  a numrer  of  curves  to  the  entire  popu- 
lation of  loading  data,  we  decided  that  it  would  be  better  to  present 
the  data  using  one  analytical  expression  for  the  range  0-10  kbar  and 
another  lor  the  entire  set  of  data  going  from  zero  to  about  33  kbar, 
Equation  (2). 

Upper  Limit  of  Utility — The  transition  from  face-centered 
cubic  (FCC)  to  body-centered  cubic  (DCC)  ytterbium,  with  its  attendant 
drastic  change  in  electrical  resistance,  begins  as  low  as  about  32  kbar 
under  dynamic  loading  conditions  (rise  time  under  500  nsec)  and  is  sub- 
stantially complete  in  less  than  2 usee.  We  have,  therefore,  concluded 
that  30  kbar  is  the  upper  limit  for  use  of  ytterbium  as  a dynamic  stress 
transducer  if  some  margin  of  safety  is  desired  by  the  experimenter. 

31  3 Divergent  Flow  Effects— At  stress  levels  around  1 kbar  and 
below  and  lor  long  times,  it  is  possible  that  changes  in  resistance  arising 
from  the  effects  of  divergence  in  the  How  field  could  be  of  sufficient 
magn  tude  to  signi f ieiantl y ailect  the  shape  of  tile  recorded  signal.  How- 
ever, no  definite  conclusions  can  be  drawn  from  our  results  because  (1) 
wo  do  not  know  how  much  tangential  strain  is  coupled  into  tho  gage  element 
in  the  field  and  (2)  the  experiments  that  were  performed  to  ascertain  tiie 
magnitude  of  strain  effects  have  shown  that  the  response  seems  to  be 
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sensitive  to  the  experimental  configuration,  suggesting  that  a gage  in 
a field  package  might  behave  differently  because  oi',  lor  example, 

Poisson  contraction  effects. 

3.1.4  Reproducibility  and  Accuracy — A variation  in  the  metallurgical 
history  of  ytterbium  foil  (annealed  vs  cold-worked)  has  a profound  effect 
on  the  magnitude  of  the  resistance  change  with  stress  at  lower  stresses 
(up  to  about  5 kbar)  and  we  have  surmised  that  the  effect  persists  at 
higher  pressures.  We  have  not  shown  any  such  effect  or  variations  in 
chemical  purity. 

W'e  believe  that  the  state  of  cold  work  of  the  foil  is  the  major 
variable  affecting  the  reproducibility  of  the  calibration  data  and  the 
accuracy  of  measurements  made  with  ytterbium  gages. 

3.1.5  Time  Effects — Xo  time  effects  of  importance  were  observed 
in  the  time  regime  of  interest  to  us  in  dynamic  measurements.  However, 
the  static  response  of  ytterbium  is  time-dependent  during  release  over 
long  periods  of  time  (minutes  to  hours)  and  should  be  used  with  caution 
for  measurements  in  static  experiments. 

3 . 2 Recommendations 

Ytterbium  shows  very  high  sensitivity,  is  well-behaved  and  sufficiently 
reproducible  (although  not  linear)  in  its  calibration,  shows  relatively 
small  sensitivity  to  noncompressi ve  stress  fields  over  most  of  its  usable 
stress  range,  is  fabricable,  and  reasonably  economical.  We  believe,  therefore, 
that  cold-rolled  ytterbium  foil  is  suitable  lor  use  as  the  active  element 
in  stress  transducers  designed  to  measure  dynamic  geological  properties 
in  situ  and  has  labiratory  applications  as  well,  such  as  making  stress 
history  measurements  in  inhomogeneous  materials.  This  section  contains 
our  specific  recommendations  and  suggestions  for  using  ytterbium  gages 
and  analyzing  experimental  data. 
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3.2.1  Ytterbium  Gages — Suggestions  concerning  the  fabrication  and 


use  of  ytterbium  foil  gages  are: 

(1)  He  sure  that  the  foil  used  is  cold-rolled  (not  annealed 

or  rolled  with  intermediate  anneals)  from  a cast  ingot--not 
from  distilled  lump  material. 

(2)  Do  not  use  ytterbium  gages  at  expected  stress  levels  above 
30  kbar  without  a back-up  system  (such  as  a manganin  or 
carbon  gage  or  a particle  velocity  gage)  in  case  the  phase 
transformation  begins  during  loading  and  makes  the  recorded 
data  unreliable. 

(3)  Do  not  expose  the  foil  to  excessive  heat  (above  about  100°C) 
for  any  length  of  time.  This  suggestion  is  made  to  avoid 
any  annealing  and  the  attendant  changes  in  sensitivity. 

(1)  Given  the  state  of  our  knowledge  about  sensitivity  to 
divergent  flow,  we  suggest  that  attempts  be  made  to 
minimize  tangential  flow  and  bending  rather  than  correcting 
for  it. 

3,2.2  Calibration  Data — We  have  not  yet  ascertained  the  best  way 
to  use  the  loading  and  release  data  presented  in  this  report,  or  if  there 
actually  is  a single  best  way  for  all  situations.  Further  experience  in 
the  systematic  application  of  the  data  would  be  necessary  before  detailed 
recommendations  for  ytterbium  users  could  be  made  in  some  type  of  handbook 
form.  However,  some  general  suggestions  about  the  use  of  the  calibration 
data  are: 

(1)  If  the  peak  stress  level  associated  with  the  event  is  10  kbar 
or  below,  use  Equation  (1)  for  the  loading  and  the  unloading  curves 
described  by  Equation  (3).  There  are  two  ways  to  do  this — graphically 
and  analytically. 

For  the  graphical  method,  plot  Equation  (1)  (shown  in  Figure  33) 

and  find  the  peak  stress  (g-  ) reached  in  the  experiment.  Then  find  the 

P 

zero-stress  residual  resistance  (AH  /R  ) in  Figure  40  and  draw  a straight 

11  o 

line  between  n and  (AR  /R  ) . This  is  the  unloading  path  from  that  peak 
P Ho 

stress.  Repeat  the  process  for  subsequent  loading  and  release  cycles. 

This  will  cause  some  error,  but  we  do  not  have  sufficient  data  to  assess 
its  magnitude. 
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The  equations  describing  the  0 to  10-kbnr  behavior  will  lend 
themselves  to  a simple  computer  code  treatment.  Essentially,  the 
computer  finds  c(t)  in  the  wave  front  using  Equation  (1)  until  it 
comes  to  a maximum.  It  knows  that 
AK 

a = A f B (— — ) (3) 

u R 

o u 

and  it  knows  Ale  ) and  B((j  ) from  the  data,  to  it  solves  Equation  (3) 

P P 

and  analyzes  the  entire  stress-time  record. 

(2)  If  the  peak  stress  is  above  10  kbar,  the  situation  is  not  as 
straightforward.  Use  Equation  (2)  for  the  loading  if  the  fit  and 
detailed  behavior  of  the  wave  around  one  kbar  and  below  are  not  of 
extreme  importance.  If  the  peak  stress  is  above  10  kbar  and  low-stress 
behavior  is  important,  we  suggest  some  combination  of  Equation  (1)  and 
some  other  function,  possibly  Equation  (2).  We  have  not  yet  found  a 
best  solution  to  this  problem  for  all  cases. 

The  unloading  data  from  stresses  above  10  kbar  should  be  used 
carefully  because  the  curves  that  represent  the  individual  data  sets 
sometimes  do  not  intersect  the  loading  curve  at  the  peak  stress  of  the 
experiment  in  which  they  were  measured.  This  often  happens  with  fitted 
data,  and  we  suggest  graphical  or  analytical  interpolation  between  the 
curves  representing  the  data. 


1 . TECHNICAL  DISCUSSION 


■1 . 1 Piezoresisti vity  and  Pic/.oresi.stanco 

The  resistivity  D ( of  a solid  is  an  intensive  property  and  is 
1 K 

representable  by  a second-rank  tensor  relating  electric  field  intensity 

E and  electric  current  density  ) , i.e., 
i k ’ 


The  resistivity  is  a function  of  the  instantaneous  state  of  strain 
existing  in  the  solid  and  of  the  strain  history.  A change  in  resistivity 

* 

caused  by  a change  in  the  state  of  strain  is  called  the  elastoresistivity 
effect.  If  we  wish  to  consider  the  change  in  resistance  to  be  caused  by 
the  application  of  stress,  we  call  the  effect  piezoresistivity.  It  has 
been  suggested  that  interpretation  of  these  effects  in  terms  of  funda- 
mental parameters  of  the  material  is  more  conveniently  formulated  with 
respect  to  the  elastoresistivity  effect,  and  measurement  problems  are 
more  easily  discussed  from  the  viewpoint  of  piezoresistivity  (Keyes,  1960). 
The  fractional  change  in  resistivity  induced  by  the  application  of 

a generalized  stress  (7  is  given  by 

mn 


ikmm  ^mn 


(6) 


where  the  77  are  the  piezoresistivity  coefficients,  which  form  a 
i Kinn 

fourtli  rank  tensor.  In  the  simplified  matrix  notation  this  can  be  written 


Valid  only  in  the  elastic  region  of  the  solid. 

It  will  also  be  shown  that  we  are  interested  only  in  principal 
directions  in  Appendix  A. 
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da 

1 


rile  resistance  of  a material  depends  not  only  on  the  resistivity, 
bul  also  on  the  dimensions.  A differential  change  in  resistance  is 
given  by  (from  differentiating  K = 0,»/A) 


cll{.  da  d(  dA 

i _i  ^ i i 

K 0 + t A 


(8) 


or 


71 . .a . 

ij  j 


(9) 


where  i corresponds  to  the  x,  y,  or  z direction,  * is  the  length 
ol  the  conductor,  A the  cross-sectional  area  normal  to  the  current 
Llow  direction,  and  11  is  the  resistance.  in  ytterbium  the  change  in 

resistivity  is  by  far  the  dominant  effect,  but  this  varies  from  material 
to  material. 


When  choosing  a material  for  the  active  element 
the  important  properties  are: 


in  a stress  transducer 


(1)  Sensitivity  The  resistance  change  with  stress  should  be 
large  enough  to  produce  a usable  signal,  and  the  dependence 
ol  the  resistance  on  stress  should  increase  (or  decrease) 
monotonically . 

(2)  Reproducibility— The  uncertainty  in  the  stress  measurements 
traceable  to  variations  in  the  gage  material  should  be 
within  acceptable  1 limits. 

(3)  Temperature  sohsi  U vi  ty-The  ratio  of  resistance  change 
induced  by  temperature  (cither  from  shock  or  joule  heating) 
to  that  arising  from  mechanical  stress  should  be  as  small 
as  possible. 

(d)  Sensitivity  to  strains  other  than  those  caused  tv  the  principal 
compressive  stress— When  a transducer  element  experiences  a 
divergent  flow  field,  there  arc  strains  induced  other  than  the 
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principal  compressive  stress.  Ic  is  important  that  the 
sensitivity  of  the  gage  material  to  stresses  other  than 
the  principal  compressive  stress  be  as  small  as  possible 
relative  to  the  compressive  stress  sensitivity.  Signals 
arising  l'rom  bending  or  tensile  stresses  can  distort  the 
gage  output  if  they  become  large  enough. 

(5)  Unloading  behavior — It  would  be  advantageous  if  the 
loading  and  unloading  behavior  of  a transducer  element 
were  the  same;  however,  unloading  behavior  can  also  be 
calibrated  experimentally,  as  we  have  done  here. 

(6)  Fabrication  and  economics — Ease  of  fabrication,  including 
ease  of  handling,  chemical  inertness,  etc.,  and  reasonable 
price  are  of  greater  or  less  importance  depending  on  the 
quantity  ol  measurements  that  are  envisioned. 

•1 . 2 Ytterbium 

Ytterbium  is  a rare  earth  element  located  between  thulium  and 

lutetium  in  the  Lanthanide  series.  Its  atomic  weight  is  173.04,  and 

o 

its  atomic  number  is  70,  At  25  C and  one  atmosphere  pressure,  its 
equilibrium  crystal  structure  is  FCC.  The  electrical  resistance  of 
ytterbium  increases  markedly  between  atmospheric  pressure  and  a maximum 
located  somewhere  between  30  and  40  kbar  (depending  on  the  state  of 
strain)  where  the  resistance  decreases  to  about  80  percent  of  the 
original  value  at  one  atmosphere.  The  maximum  resistance  change  under 
hydrostatic  pressure  has  been  reported  to  be  as  much  as  23  times  the 
value  of  the  resistance  at  one  atmosphere  (Stromberg  and  Stephens,  1964) 
and  as  little  as  6 times  (Stager  and  Drickamer,  1963).  This  phase 
transition  is  especially  interesting  because  the  high  pressure  modifi- 
cation is  BCC , which  is  usually  more  open  than  the  low  pressure  FCC 
structure,  but  the  phase  change  does  involve  an  overall  volume  decrease 
(Stephens,  1964). 
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The  phase  change  has  been  discussed  by  Hall  and  Merrill  (1963), 
lhill  et  al.  (1963),  Souers  and  Joura  (1963),  and  Stager  and  Drickamer 
(1963).  All  agreed  that  the  FCC  - BCC  transformation  was  accompanied 
by  a decrease  in  the  size  of  the  individual  ytterbium  atoms,  and  the 
necessary  shift  in  electronic  structure  was  the  cause  of  the  decrease 
in  resistance  and  return  to  metallic  conduction  behavior.  However, 
the  kinetics  of  I he  conductivity  change  do  not  fit  with  an  atom  shrinkage 
model.  Jhis  is  c iseussed  in  more  detail  in  Section  4,7,6. 

Other  polymorphic  changes  that  have  been  noted  are  a FCC  — BCC 


transi ormation  at  one  atmosphere  and  approximately  798°C  (Daane,  1961),  and 


a shear  transformation  from  a FCC  to  hexagonal  close-packed  (HCP)  structure 

o 


occurring  at  one  atmosphere  and  approximately  0°C  (Kayser  and  Soderquist, 
1969,  Kayser,  1970;  Tanurna  et  al.,  1970;  Bucher  et  al , , 1970).  Results 
showing  other  temperature-induced  effects  have  been  presented  by  Hurd 
and  Alderson  (1973),  The  large  difference  between  the  extremes  of  the 
peak  resistance  values  measured  by  various  investigators  might  have  arisen 
from  impurity  eifects  at  high  pressures,  from  the  effect  of  nonhydrostatic 
components  of  stress,  or  from  differences  in  sample  preparation  techniques. 
As  the  investigators  who  have  studied  this  proolcm  have  improved  their 
techniques,  they  have  obtained  fractional  resistance  changes  at  the 
transition  stress  which  are  in  the  range  of  6 to  8. 

It  is  not  surprising  that  a substance  that  displays  such  complex 
behavior  has  attracted  considerable  attention  from  both  experimentalists 
and  theoreticians.  The  large  dependence  of  electrical  resistance  on 
stress  has  also  aroused  a practical  interest  in  the  properties  of 
ytterbium  because  of  it?  attractiveness  for  use  as  a quantitative  dynamic 
stress  sensor,  especially  in  the  low  stress  (1-20  kbar)  region.  Develop- 
ment and  calibration  programs  have  been  recently  carried  out  at  SRI 
(Keough  1970;  Williams,  1970;  Ginsberg,  1971),  and  ytterbium  gages  have 
been  used  successfully  in  laboratory  experiments  (Petersen  and  Erlich, 

1972;  Erlich  et  al.,  1971)  as  well  as  in  the  field. 
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^ ^ Pieparntion  and  Characteri zat  ion  of  Specimens 

1 t terbium  metal  is  obtained  by  mixing  its  oxide  with  a reducing 
agent  such  as  calcium,  silicon,  aluminum,  or  magnesium  or  a combination 
oi  these  metals)  and  heating  the  mixture  under  vacuum.  The  oxide  reacts 
with  the  reducing  agent  and  the  resulting  ytterbium  vapor  is  condensed 
on  a cold  collector  placed  above  the  crucible.  The  metal  obtained  from 
this  process  is  often  purified  further  by  redistillation  under  vacuum, 
and  the  final  product  is  called  distilled  lump.  This  lump  can  be  rolled 

directly  into  foil,  or  it  can  be  remelted,  cast  into  ingot,  and  then 
rolled  into  foil. 

The  specimens  used  in  our  dynamic  calibration  experiments  were 
tour-terminal  grids  (see  Figure  1)  fabricated  from  as-rolled  (from  ingot) 
ytterbium  toil  of  nominal  0.005-cm  (2-mil)  thickness.  Most  of  the 
experiments  were  done  on  foil  purchased  from  Research  Chemicals,  Inc., 
Phoenix,  Arizona.  For  comparison  purposes,  some  work  was  done  with  foil 
obtaind  from  American  Rare  Earth  and  Foil  Co.,  Tempe,  Arizona.  Both 
suppli  claimed  99.9%  purity.  A spectrographic  analyses  of  five  sample* 
for  the  Research  Chemicals  (RC)  material  is  shown  in  Table  1.  The  purity 
of  the  material  from  American  Rare  Earth  and  Foil  (AREF)  was  also  quoted 
to  be  99.9%.  An  analysis  is  included  in  Table  1. 

Each  grid  was  made  from  a 2-cm  square  cut  from  the  as-received  foil. 
The  thickness  of  all  squares  was  measured,  they  were  cleaned  carefully 
with  a small  glass  fiber  brush,  and  the  grid  photoetched  on  an  epoxy- 
fiberglass or  strippable  mylar  substrate  using  a mixture  of  ethanol 
and  hydrochloric  acid.  The  cover  of  fiberglass  cloth  was  added  and 
then  also  impregnated  with  epoxy,  making  the  package  shown  in  Figure  2. 
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Table  1 

SPECTROGRAPH I C ANALYSIS  OF  YTTERBIUM  FOIL 


Sample  No. 


Element 

RC-1 

RC-2 

RC-3 

RC-4 

RC-5 

AREF-1 

Ytterbi urn 

Major 

Major 

Major 

Major 

Major 

Major 

Thulium 

<0 . 005 

<-0 . 005 

<0 . 005 

<0.005 

<0.005 

ND 

Erbium 

0.006 

0.02 

0.005 

<0.005 

0.005 

ND 

Lutetium 

<0 . 005 

<0.005 

■'0.005 

>-0 . 005 

<0.005 

ND 

Dysprosium 

<0.005 

<0.005 

<0.005 

<0.005 

<0.005 

NA 

Yttrium 

<D.005 

<0.005 

<-0.005 

<-0.005 

<-0.005 

NA 

Silicon 

0.01 

0.01 

0.C1 

0.01 

<0.01 

ND 

I ron 

0 . 01 

0 . 01 

0.02 

0.03 

0.02 

ND 

Tantalum 

ND 

ND 

ND 

ND 

ND 

* 

<0.5 

Magnesium 

0.01 

0.01 

0.01 

0.01 

0.01 

-0.001 

Calcium 

-o.oi 

<-0.01 

<0.01 

<-o.oi 

<-o.oi 

<0.001 

Aluminum 

<0.01 

0.02 

<-o.oi 

0.01 

0.02 

<-0.001 

Nickel 

NA 

NA 

NA 

NA 

NA 

ND 

Copper 

NA 

NA 

NA 

NA 

NA 

ND 

Manganese 

NA 

NA 

NA 

NA 

NA 

■-0.001 

NU  = not  detected 
NA  = not  analyzed  for 
< = less  than 


No  standards  for  Ta  content  less  than  this  were  available. 
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SHOCK  DIRECTION 


FIGURE  2 YTTERBIUM-FIBERGLASS-EPOXY  PACKAGE  (EXPLODED  VIEW) 

4-3*1  Resistivity  Measurements — We  determined  the  resistivity  of 
a number  of  samples  of  both  the  RC  and  AREF  foils  by  carefully  measuring 
the  line  lengths  and  widths  of  individual  grids  (the  thickness  had  been 
measured  previously)  and  then  measuring  the  resistance  with  a Cimron  Model 
6853  digital  multimeter.  The  results  are  contained  in  Table  2. 


Table  2 

RESISTIVITY  OF  0.005-cm  (2-mil)  YTTERBIUM  FOIL 


Foil  Type 


RC  (Lot  109-C,  as  rolled) 
AREF  (as-rolled) 

AREF  (annealed) 


Resistivity  (p) 
p,Q-cm  

27  ± 3 

39  ± 5 

19  ± 5 
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The  value  of  resistivity  usually  quoted  for  ytterbium  is  around 
30  uQ-cm.  Higher  values  have  been  reported  (Keough,  1970;  Williams,  1970), 
but  the  value  of  p we  obtained  for  the  annealed  AREF  foil  is  lower  than 
any  ve  have  seen  reported.  The  major  uncertainty  is  in  the  thickness 
measurement.  Although  we  know  that  the  electrical  properties  of  ytterbium 
are  strongly  dependent  on  its  metallurgical  history,  its  resistance  should 
also  be  very  sensitive  to  its  purity,  so  we  cannot  unequivocally  account 
for  the  differences  among  the  three  readings,  except  that  the  low  value 
foi  the  annealed  foil  would  seem  to  be  qualitatively  correct. 

1 • 3 • 2 Measurements  of  Temperature  Coefficient  or  Resistance— We 

measured  the  temperature  coefficient  of  resistance  of  four-terminal 

ytterbium  gages  packaged  in  thin  fiberglass  by  immersing  them  in  a 

silicone  oil  bath  and  measuring  the  resistance  change  of  the  gage  as  a 

function  of  temperature.  The  temperature  of  the  oil  bath  was  monitored 

with  mercury  thermometers;  the  bath  was  insulated  and  stirred  to  minimize 

temperature  gradients  within  the  bath.  The  resistance  of  the  gages  was 

measured  with  a Cimron  Model  6853  digital  multimeter  having  a specified 

accuracy  of  0.02%  in  the  range  of  interest.  Our  own  cross-checks  against 

a recently  calibrated  Kelvin  double  bridge  confirmed  the  accuracy. 

The  power  dissipated  in  the  gage  by  the  ohm  meter  current  was  only  ^0.01 

watt,  and  the  temperature  of  the  gage  did  not  rise  (evidenced  by  an  absence 

of  resistance  drift)  during  the  course  of  a measurement. 

The  RC  foil  had  a temperature  coefficient  of  approximately  1.35  x lo”3 

VQ  C,  between  30  and  100  C,  and  the  AREF  material  had  a lower  coefficient 

“3  * 

approximately  0.51  x 10  Q/Q  C.  The  annealed  material,  also  AREF,  had 

a higher  coefficient  than  either  of  the  other  two  as-rolled  foil  samples, 

this  being  2.65  x 10  Q/Q  C.  Keough  (1970)  and  Ginsberg  (1971)  also 

found  that  higher  resistivity  material  tends  to  have  a lower  temperature 

coefficient  of  resistance,  but  we  are  not  able  to  present  an  explanation 

for  this  observation  based  on  our  experimental  data. 
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l.*5.3  Hard  ness --The  hardness  of  representative  pieces  of  foil  was 
measured  with  a Kentron  microhardness  testing  machine,  using  a diamond 
py 1 amid  indenter  and  a 10— gram  load.  Six  indentations  were  made  in  each 
2 x 2 cm  square . We  were  hoping  that  these  tests  would  be  a good  foil 
quality  control  test  because  hardness  is  a good  measure  of  the  amount  of 
cold  work  a metal  has  ex per.  enced.  We  found,  however,  that  the  measured 
hardness  of  the  annealed  material  did  not  differ  from  the  measured 
hardness  of  the  as-rolled  material,  which  is  very  unlikely.  Our  average 
measured  DPN  hardness  for  the  two  lots  of  as-rolled  foil  was  23  ± 2. 

Although  performing  hardness  measurements  on  thin  foils  is 
difficult,  we  at  least  expected  that  the  test  could  distinguish  between 
annealed  and  unannealed  foil.  Our  inability  to  do  so  was  ascribed  to  the 
presence  of  a surface  oxide  film  on  the  foils.  Mechanical  removal  of  the 
film  hardens  the  surface,  thereby  making  hardness  measurements  meaningless, 
and  wo  were  not  able  to  find  a good  chemical  means  of  removal.  It  would 
be  very  useful  if  some  means  to  do  these  tests  accurately  could  be  devised 
because  it  represents  a possible  way  to  separate  potentially  bad  gage 
material  from  acceptable  material. 

Held  Strength  in  Tension — The  mechanical  response  of  a gage 
to  a divergent  flow  field  is  dependent  on  its  mechanical  properties  in 
tension  normal  to  the  direction  of  the  principal  stress.  We  decided, 
therefore,  that  it  would  be  useful  to  ascertain  the  yield  stress  of  our 
ytterbium  foil  gage  stock  in  simple  tension.  This  would  also  help  us 
interpret  our  tensile  strain  piczorcs: stan  experiments,  whi cli  arc 
described  later. 

Four  small  tensile  specimens  were  prepared  from  the  RC  as- rolled 

O.OOo— cm  (2— mil)  ytterbium  loil  and  tested  on  an  Instron  mechanical 

testing  machine  at  a strain  rate  of  0 . 0251-cm/inin . Three  of  the  specimens 

were  made  so  that  the  long  dimension  was  the  direction  of  the  rolling 

texture  in  the  foil;  the  other  was  made  with  the  rolling  texture  normal 
* 

Kichle  Division  of  Ametck,  Inc.,  East  Moline,  Illinois. 
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to  the  direction  of  applied  strain.  Our  best  estimate  of  the  yield 

strength  of  the  foil  used  as  gage  stock,  which  is  believed  lo  be 

representative  of  our  field  gage  stock,  was  0.50  ± 0.06  kbar.  The 

corresponding  yield  strain,  assuming  a modulus  of  187  kbar,  is  0.0027. 

This  figure  was  obtained  by  estimating  the  first  deviation  from  linearity 

* 

in  the  measured  stress-strain  data.  No  difference  was  seen  in  yield 
stress  between  the  specimens  pulled  with  the  rolling  texture  parallel 
and  normal  to  the  direction  of  applied  strain,  but  this  could  very  well 
be  attributed  to  rather  large  experimental  uncertainty. 

'1,3.5  Static  Compressive  Mechanical  Properties — Since  the  change  in 
slope  in  the  stress-resistance  curve  was  considered  a yield-point  effect 
(Ginsberg,  1971) , it  was  thought  useful  to  make  some  static  measurements 
and  obtain  an  approximate  lower  bound  value  of  the  static  compressive 
yield  point  of  ytterbium.  To  do  this,  a number  of  small  cylinders  with 
varying  diameter  to  height  (D/ll)  ratios  (Dietei  1961)  were  machined 
out  of  99.9rc  pure  ytterbium  ingot  that  was  on  hand  and  stress-strain 
tests  were  performed  on  these  cylinders,  using  an  Instron  testing  machine 
in  the  compressive  mode  of  operation.  Although  the  experiments  are  straight- 
forward, the  interpretation  of  the  results  is  not  simple.  Frictional 
effects  cause  the  pressure  distribution  to  vary  from  the  edge  to  the 
center  of  the  cylinder,  and  the  elfects  of  porosity,  which  were  undoubtedly 
present,  tend  to  lower  the  measured  vield  point.  However,  these  experi- 
ments, the  results  of  which  are  contained  in  Table  3,  have  given  us  a 
good  qualitative  estimate  of  the  yield  point  in  compression  so  that  we 
can  more  easily  discuss  the  piezoresistance  of  ytterbium  where  yield 
point  effects  could  be  important.  The  extrapolated  (D/H  = 0)  value  of 
the  compressive  yield  point  is  about  0.27  kbar,  which  we  consider  a 
lower  bound  because  of  the  porosity  and  the  fact  the  ingot  is  not  work- 
hardened  . 

* 

Engineering  stress  and  strain, 

j 

This  ingot  was  used  by  Keough  (1970)  as  a source  of  material  for 

vapor-deposited  gages. 
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Table  3 


COMPRESSIVE  YIELD  STRENGTH  OF  YTTERBIUM  INGOT 


Specimen 
No . 

Height 

cm 

Diameter 

cm 

D/ II 

C (measured) 
kb  nr 

] 

0.635 

0.917 

1.50 

o 

o 

-H 

Q 

2 

0.635 

0.792 

1.25 

0.36  ± 0.07 

3 

0.635 

0.635 

1.00 

0.30  ± 0.07 

4 

0.635 

0.321 

0.50 

0.30  ± 0.07 

•'.3.6  X 

-Ray  Crystallography — Samples 

ol  ytterbium  foil  were  examined 

by  standard 

X-ray  diffraction 

techniques , 

using  nickel-filtered  copper 

radiation  and 

a General  Electric  XRD-3  diffractometer. 

The  most  intense 

rom 


peaks  or  each  diffraction  pattern  could  be  indexed  as  reflection  fr< 

ICC  ytterbium  hating  a lattice  constant  of  5.486  A.  Most  of  the  patterns 
contained  additional  reflections,  which  could  be  indexed  as  the  strongest 
reflections  of  IiCP  ytterbium  having  lattice  constants  c = 6.328  A and 

3.884  A.  The  relative  intensities  of  IICP  and  FCC  reflections  provided 
a quantitative  measure  of  the  proportion  of  IICP  phase  in  the  ytterbium 
samples  that  were  judged  to  consist  of  randomly  oriented  fine  grains, 
on  the  basis  of  the  relative  intensities  of  FCC  diffraction  peaks.  These 
randomly  oriented  samples  contained  about  10%  IICP  and  the  proportion  of 
HCP  remained  constant  for  samples  subjected  to  plastic  deformation  and 
to  temperature  cycling  over  the  range  of  -196°C  to  150°C.  A sample  of 
shock-loaded  ytterbium  foil  was  also  examined,  and  the  proportion  of  IICP 
was  essentially  the  same  as  in  the  original  material. 

In  the  flat  sample  diffractometer  geometry  that  was  used,  the 
diffraction  pattern  originates  from  those  crystallographic  planes  that 
are  oriented  parallel  to  the  sample  plane,  i.e.,  parallel  to  the  surface 
of  the  foil.  Preferred  orientation  of  the  fine  individual  grains  that 
make  up  the  sample  may  he  inferred  by  comparison  of  observed  relative 
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diffraction  peak  intensities  with  the  relative  intensities  calculated  for 
nn  assemblage  of  randomly  oriented  grains.  We  noted  large  differences  in  the 
extent  of  preferred  orientation  among  samples  from  different  lots  of  foil. 

Ihe  criterion  for  preferred  orientation  was  the  intensity  ratio  of  (220) 
reflection  to  (111)  reflection  for  the  sample.  For  randomly  oriented 
grains  the  calculated  intensity  ratio  I (220)/I (111)  is  0.4.  One  group 
of  foil  samples,  having  I (220)/I (111)  in  the  range  of  0.3  to  0.6,  was 
classed  as  isotropic,  within  the  limits  of  the  measurement  technique. 

Another  group  of  foil  samples,  having  I (220)/I  (in ) in  the  range  of 
1.2  to  3.0  was  classed  as  having  moderate  preferred  orientation.  A 
third  group  of  samples,  for  which  I (220)/I (111)  was  in  the  range  of 
6.0  to  8.0,  was  classed  as  showing  very  strong  preferred  orientation. 

The  type  of  preferred  orientation  noted,  an  alignment  of  (110)  planes 

parallel  to  the  plane  of  the  foil,  is  common  to  all  heavily  rolled  FCC 
metals . 

HCP  reflections  were  observed  in  the  diffraction  patterns  of  the 
ytterbium  samples  that  showed  strong  preferred  orientation,  but  the 
preferred  orientation  interfered  with  quantitative  determination  of  the 
proportion  of  HCP.  Qualitatively,  we  estimate  that  the  HCP  content 
does  not  differ  significantly  from  the  10%  observed  in  the  isotropic 
ytterbium  samples.  All  of  the  samples  of  commercially  obtained  foil  and 
bulk  ytterbium  that  we  have  examined  appear  to  contain  about  10%  of  the 
HCP  phase.  The  only  samples  that  contained  no  detectable  HCP  phase  were 
films  obtained  by  vacuum  vapor  deposition  on  kapton  substrates.  These 
films  were  from  the  same  batch  as  the  samples  studied  by  Spataro  (1972). 

1 - 4 Dynamic  Calibration  Experiments 

In  the  overall  program  the  highest  priority  was  given  to  obtaining 
dynamic  loading  and  release  calibration  data  that  would  be  immediately 
applicable  to  the  analysis  of  field  data.  Less  emphasis  was  placed  on 
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the  other  parts  of  the  program,  although  they  were  important  to  our 
understanding  of  the  behavior  of  ytterbium  and  the  field  gage.  We  chose 
this  method  because  we  felt  that  a good  data  base  was  needed  before  other 
effects  ol  unknown  magnitude  could  be  sought. 

All  of  the  specimens  used  in  the  dynamic  calibration  shots  were 
encapsulated  in  epoxy-impregnated  fiberglass  packages— smaller  versions 
of  the  field  transducer  described  by  Grady  et  al.  (1973)  and  Smith  et  al . 
(1972)  and  shown  in  Figure  3.  This  was  done  because  Keough  (i970)  suggested 
that  piezorcsistant  stress  transducer  response  can  be  affected  by  the 
material  surrounding  it  and,  in  fact,  by  the  target  material  itself.  We 
did  not  wish  to  test  this  hypothesis  in  detail  because  of  the  effort 
involved,  hut  chose  rather  to  encapsulate  the  specimens  in  the  same 
immediate  environment  they  would  see  if  they  were  in  the  field,  i.e,, 
epoxy-fiberglass.  This  procedure  reduces  some  of  the  uncertainty  in 
transferring  laboratory  data  to  the  field,  but  there  was  a price  to  be 
paid  in  terms  of  response  time  and  disturbances  to  the  main  wave  arising 
from  the  presence  of  the  gage  package.  This  problem  will  be  treated  in 
more  detail  later  when  individual  shots  are  discussed.  In  the  next  sections 
we  will  describe  our  experimental  facilities  and  the  four  basic  shot 
designs  we  used  to  obtain  the  dynamic  calibration  data. 

1,1,1  jjgPcrircental  Facilities— The  dynamic  calibration  experiments 
were  performed  on  our  10.16-cm  and  6.35-cm  bore  light-gas  guns.  Helium 
was  used  as  the  driving  gas.  The  velocity  was  measured  with  pins  which 
triggered  and  then  stopped  a pair  of  digital  time  counters,  and  another 
independent  set  of  pins  which  triggered  and  produced  time  marks  on  a 
Cord in  oscilloscope.  Alignment  of  the  target  was  accomplished  by  a 
device  which  lines  up  the  target  mount  at  the  end  of  the  barrel  with  the 
gun  barrel.  The  targets  were  made  flush  with  their  respective  mounting 
rings,  and  were  installed  on  the  prealigned  target  mount.  All  recording 
was  done  with  oscilloscopes,  and  a Pulsar  Model  301  three-channel  constant 
current  power  supply  was  used  to  power  the  gages. 
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FIGURE  3 YTTERBIUM  STRESS  GAGE— RESISTANCE  GRID  ENCAPSULATED  IN  AN 
EPOXY-FIBERGLASS  MATRIX.  (UPPER)  IN  SITU  PACKAGE,  (LOWER) 

PACKAGE  FOR  GROUTING  IN  DRILL  HOLES 

A circuit  for  monitoring  the  pressure-induced  resistance  change 

in  the  transducer  is  shown  in  Figure  1.  K is  the  gage  resistance  and 

G 

R is  a resistance  in  series  with  the  recording  transmission  cable. 

b 

The  cable  is  terminated  at  the  oscilloscope  with  R , a resistance  equal 
to  the  characteristic  impedance  of  the  cable.  I.  the  constant  current, 
is  provided  by  the  gage  power  supply.  To  provide  proper  termination 
at  both  ends  of  the  cable,  the  usual  procedure  has  been  to  select  R 

S 

such  that  the  sum  of  R and  R is  close  to  the  characteristic  impedance 

S G 

of  the  cable  (Keough,  1967).  The  assumption  of  constant  current  in  the 
gage  is  not  satisfied,  however,  since  a shunting  current  ( I in  the  figure) 
occurs.  This  creates  an  error  in  the  gage  factor  and  a correction 
including  the  effect  of  the  shunt  current  must  be  performed.  When 
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FIGURE  4 CIRCUIT  FOR  MONITORING  TRANSIENT  RESISTANCE 

CHANGES  IN  YTTERBIUM  FOILS 


ytterbium  is  used  as  the  transducer  material  the  correction  term  is 
considerable  and  the  complication  can  be  lessened  by  adjusting  R such 
that  the  constant  current  assumption  is  satisfied.  Single-ended  .able 
termination  is  adequate.  In  the  present  work  values  of  equal  to  500 
and  1000  i.  were  used  at  higher  stress  levels.  This  precaution  reduced 

S,Umtl"B  CUr,'ent  t0  **-  »/**  •*  the  total  current  1»  individual 

experiments.  The  voltage-time  records  were  digitized  directly  Iron,  the 
oscillograms  with  a Whittaker  Telocordex  counter.  The  output  fed 
into  a keypunch  and  the  cards  were  used  to  plot  the  data  and  tabulate 
the  fractional  voltage  change  as  a function  of  time.  Conversion  of 
IV  i;  to  AH/!io  was  done  using  Wilkinson's  analysis  (in  Keough,  1968)  of 
the  recording  circuit  which  is  discussed  in  more  detail  in  Append!* 

■‘■4'2  jltet-Stress  Step-Release  Experl ments-The  step-release 
technique  (Keough,  1970)  produces  one  peak  stress  calibration  point  and 
a series  of  points  on  the  unloading  path  originating  a,  the  peak  stress 
obtained  during  the  experiment.  When  a thin,  high-impedance  flying  plote 
impacts  a thicker  target  of  lower  shock  impedance,  a location  in  the 
target  near  the  impact  surface  will  firs,  ioad  „p  to  so>ie  peak 

and  then  experience  a series  of  transmitted  release  waves  originating 
at  the  back  of  the  projectile  head.  If  both  the  projectile  material 
and  the  target  behave  elastically  (linear  or  nonlinear)  in  the  stress 
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range  ol  interest,  a series  of  stress  states  on  one  unloading  path  can 
be  determined  i rom  a stress-particle  velocity  diagram  of  the  two  materials. 
Figure  6 shows  a schematic  of  typical  experiment  design,  and  Figure  7 is  a 
schematic  x-t  diagram  which  depicts  the  expected  wave  shape.  Figure  8 is 
an  impedance  match  diagram  showing  the  stress  states  reached  in  a step- 
release  experiment. 

The  targets  were  constructed  from  two  pieces  of  10*15-c.m  diameter 
GE  Type  151  fused  quartz.  The  dynamic  properties  of  this  material  have  been 
studied  in  detail  by  Barker  and  hollenbach  (1970).  The  front  pieces 
(impact  surfaces)  were  0 . 2 cm  thick  and  the  backs  approximately  1.9  cm 
thick.  The  gage  leads  were  brought  out  through  holes  drilled  in  the  back 
piece.  In  order  to  improve  the  probability  of  lead  survival,  dental 
amalgam  was  packed  into  the  holes  after  the  copper  leads  were  brought  out. 
This  lead  arrangement  was  very  successful,  and  recording  times  of  up  to 
10  usec  were  achieved  at  stress  levels  of  up  to  20  kbar  and  m-re.  The 
gage  packages  were  glued  to  the  back  piece  with  Homolite  epoxy.  More 
epoxy  was  then  poured  around  the  packages  and  the  front  plate  was  pressed 
on,  giving  a total  gage  plane  glue  li  ie  and  package  thickness  of  between 
0.020  and  0,025  cm.  Figure  9 is  a plot  rf  a digitized  oscillogram  taken 
during  Shot  S19,  the  records  from  which  are  shown  in  Figure  10.  The  steps 
are  evident,  which  allowed  us  to  obtain  good  release  data,  but  it  can  be 
seen  that  there  were  some  transient  disturbances  to  the  flow  caused  by 
the  presence  of  the  low-impedanc.e  gage  packages  and  glue  line.  There 
was  also  some  apparent  dispersion,  which  was  caused  by  a small  amount  of 
shock  tilt  in  conjunction  with  the  presence  of  the  thick  glue  line  and 
gage  package. 

* Or,  if  (as  shown  in  Figure  5)  the  release  paths  from  all  stress  states 
up  to  and  including  the  highest  peak  stress  in  the  experimental  series 
lie  on  the  Hugoniot  and  are  coincident,  as  in  Plexiglas  up  to  about  6 kbar. 
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STRESS  — kbar 


FIGURE  5 UNLOADING  STRESS-STRAIN  DATA  FOR  PMMA  (BARKER  AND  HOLLENBACH, 


PROJECTILE:  LUCALOX  (HOT-PRESSED  ALUMINUM  OXIDE)  HEAD,  0.5  cm  THICK 
6061-T6  ALUMINUM  BODY 
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GA-7511-12A 


FIGURE  8 RELEASE  STRESS  STATES,  IMPACT  OF 
HIGH  SHOCK  IMPEDANCE  MATERIAL 
ON  LOWER  IMPEDANCE  TARGET 


MP  1797-13 

FIGURE  10  OSCILLOSCOPE  RECORDS  FROM  SHOT  S19 
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I'liLf  power  supply  and  slow-sweeping  oscilloscopes  were  triggered  with 
radial  pins  placed  2 to  5 mm  from  the  impact  surface,  the  arrival  of  the 
projectile  closing  a circuit  between  a pin  and  ground.  Faster  sweeping 
oscilloscopes  were  triggered  with  copper  foil  switches  placed  on  the 
impact  surface. 

Two  types  of  projectiles  were  used,  depending  on  the  velocity  required 
to  reach  the  proper  stress  level.  Both  were  fabricated  from  6061-T6  aluminum 
alloy  and  had  10.15  cm  diameters,  but  the  projectile  used  for  the  lower 
stress  (8-15  kbar)  experiments  weighed  approximately  15  kg  and  was  75  cm 
long,  while  the  projectile  used  above  this  stress  range  weighed  approxi- 
mately 2 kg  and  was  18  cm  long.  The  projectile  heads  were  0.5-cm  thick 
Lucalox  plates  approximately  9.9  cm  in  diameter.  These  plates  were  ground 
flat  and  parallel  to  within  0.0013  cm  and  were  glued  to  the  projectile  with 
epoxy  cement,  Ihere  was  a 0.635-cm  hollow  behind  each  head  so  that  the 
x elease  could  originate  at  the  rear  free  surface  of  the  projectile  head. 

1 liese  shots  were  designed  for  a 4.5-gsec  one-dimensional  strain 
state  in  the  projectile  head  over  the  gage  package  area.  Actual  time  of 
one-dimensional  flow  in  the  target  was  greater  than  this,  so  that  data 
obtained  at  times  greater  than  4.5  ^sec  were  accepted  as  valid. 

Low-Stress  Step-Release  Experiments — Figure  11  is  a schematic 
diagram  of  the  target  and  projectile  used  for  a series  of  step-release 
| shots  in  the  0 . 75-to-6-kbar  range.  The  targets  were  made  from  two 

pieces  of  7 . 62-cm  diameter  Type  II  UVA  Plexiglas/  The  dynamic  mechanical 
i;  properties  of  this  material  have  also  been  studied  by  Barker  and  Hollenbach 

(1970).  The  front  pieces  (impact  surfaces)  were  0.2  cm  thick,  and  the 
backs  were  0.312  cm  thick.  The  gage  leads  were  brought  out  the  sides  of 
the  targets.  Recording  times  of  10  p,sec  and  more  were  routinely  obtained 
vith  this  configuration.  The  gage  packages  were  glued  into  a slot  milled 
into  the  rear  section,  and  then  the  target  was  assembled  with  Homolite  epoxy. 

General  Llectric  Co.  brand  of  hot-pressed,  high  density  aluminum  oxide, 

1 Rohm  and  Haas,  In  . 
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O-RING  GROOVES 
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6.32  cm 


GLASS 

HEAD 


PROJECTILE:  MIRROR  OR  SODA-LIME  GLASS  HEADS.  EITHER  0.318  cm.  0.475  cm.  OR  0.635  cm  THICK 
6061-T6  ALUMINUM  BODY 


— | |— — 0.15  cm 


1.91  cm 

TARGET:  TYPE  II  UVA  PLEXIGLAS.  DIAMETER:  7.62  cm 
Fiberglass-epoxy  gage  package  is  approximately  2 cm  square  and 
0.025  cm  thick.  Gage  grids  are  0.635  cm  square  with  ten  elements, 
each  0.005  cm  thick  and  0.05  cm  wide. 

MA-1 797-7 

FIGURE  11  SCHEMATIC  DIAGRAM  OF  GLASS-PLEXIGLAS  STEP-RELEASE  EXPERIMENTS 
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Hie  total  thickness  of  the  gage  packages  and  glue  was  about  0.025  cm. 

The  oscilloscopes  and  power  supply  were  triggered  in  the  same  way  as  the 
other  step-release  shots. 

Figure  12  is  a plot  of  the  data  obtained  during  Shot  Sid,  and 
Figure  IS  shows  the  actual  records.  Because  the  impedance  of  the  epoxy 
and  the  epoxy-impregnated  fiberglass  is  closer  to  that  of  Plexiglas  than 
it  is  to  fused  quartz,  the  transient  disturbances  are  of  much  smaller 
magnitude  and  do  not  interfere  in  any  important  way  with  the  wave  shape. 

The  amount  of  dispersion  is  very  small,  and  the  data  from  this  series 
are  generally  of  very  good  quality  with  easily  discernible  steps. 

Two  types  of  projectiles  were  again  used,  the  choice  depending  on 
the  velocity  required  in  the  experiment.  For  low  stress  shots  (2  kbar 
and  below)  we  used  30.5-cm  long,  6.32-cm  diameter  hollow-bodied  projectiles 
weighing  approximately  2.2  kg.  To  achieve  stress  levels  between  2 and 
about  G kbar,  we  used  a 15.25-cm  long,  6.32-cm  diameter  hollow-bodied 
projectile  which  weighs  approximately  0.6  kg.  The  projectile  heads 
were  glass  disks  either  0.318,  0.476,  or  0.G35  cm  thick  with  5.75-cm 
diameters.  The  disks  wore  glued  into  a recess  in  the  front  of  the 
projectiles  with  epoxy  cement.  A 0.158-cm  vacuum-relieved  hollow  area 
was  left  behind  each  disk.  The  density  and  sound  velocity  of  each  disk 
were  measured  separately  so  the  peak  stress  and  the  steps  could  be 
accurately  calculated. 

I he  maximum  stress  at  which  this  shot  configuration  could  be  used 
was  approximately  G kbar.  There  are  no  rate  effects  below  this  stress 
level,  and  the  release  paths  for  Plexiglas  obtained  by  Barker  and  Ilollenbach 
(1970)  no  longer  coincide  with  the  Uugonir.t  of  the  material  nor  do  they 
coincide  with  each  other  above  this  limit.  Those  conditions  of  no  rate 
effects  and  coincidence  of  Ilugoniot  and  release  paths  w'ere  required  for 
calculation  of  the  release  states  in  the  calibration  experiments.  The 
glass  was  assumed  to  be  linearly  elastic  in  this  low  stress  region 
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FIGURE  13 


OSCILLOSCOPE  RECORDS  FROM  SHOT  S14 
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<SU!""an  Ct  a1,  1909  = KoSOnborK  ilnci  fiinsborK|  1972)  and  the  two  types  used 
h:K'  Khock  imP«-Mlance.s  of  11.7  ± 0.1  and  11.5  ± 0.1  g . mm/cm3  . soc 

UV  expected  that  the  minimum  one-dimensional  strain  recording  time, 
calculated  by  finding  the  time  of  arrival  of  edge  rarefactions  at  the 
gage  area  using  the  glass  projectile  head  shook  velocity  rather  than 
the  Plexiglas  target  velocity,  was  about  1 gscc.  Actual  usable  recording 
time  was  considerably  greater  than  this. 

,'1’1  ■1~-Hblc’~‘C;h()C'k  fxpe  ri  men  ts— Figure  11  is  a schematic  diagram 

of  the  target  and  projectile  used  for  a series  of  low  stress  experiments 
in  which  we  obtained  two  loading  points  and  data  on  the  zero-stress 
residual  resistance  in  each  experiment.  Figure  15  is  a schematic  diagram 
ol  the  wave  interactions,  and  Figure  10  shows  how  the  second  wave  stiess 
level  is  obtained. 

The  targets  were  again  made  from  two  pieces  of  7.02-cm  diameter 
I.vpo  II  LVA  Plexiglas.  The  front  piece  was  0.15  cm  thick  and  the  back 
°’:,17  Cm  thick‘  C°PP°r  Metrical  leads  were  brought  out  the  sides  of 
the  targets.  Recording  times  of  up  to  10  usee  were  obtained  at  the 
lowest  stress  levels.  The  targets  were  assembled  in  the  same  way  as  the 
i°w  StrOSS  •steP-rGloase  shots  described  earlier,  and  power  supply  and 
oscilloscope  triggering  was  also  accomplished  by  the  same  methods. 

The  projectiles  were  6061-T6  aluminum,  50.5  cm  long  and  6.52  cm  in 
diameter,  and  weighed  about  2.2  kg.  The  projectile  heads  were  1.06-cm 
thick  disks  of  Plexiglas  glued  to  the  projectile  with  epoxy  cement. 

There  was  no  hollow  space  behind  the  head,  so  stress  release  originated 
at  the  rear  free  surface  of  the  target. 

17  a plot  of  the  da  a obtoined  during  Shot  S7  ond  Figure  18 
contains  the  octuai  records.  The  small  spike,  occurring  between  5 and 

6 ■— * 0ml  l,0tWeon  7 "n<:  » “*«>  from  flow  disturbances  caused  by 

the  gage  package,  but  the  trace  returns  to  t.,o  e.,uilibri„„  level  very 

duickly.  We  consider  the  value  of  at  the  first  peak  and  the  residual 

resistance  at  aero  stress  (SK^IO  to  be  more  reliable  data  than  the 
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PROJECTlLt:  TYPE  II  UVA  PLEXIGLAS  HEAD,  1.06  cm  THICK 
6061-T6  ALUMINUM  BODY 


TARGET:  TYPE  II  UVA  PLEXIGLAS.  DIAMETER:  7.62  cm 
Fiberglass-epoxy  gage  package  is  approximaiely  2 cm  square  and 
0.025  cm  thick.  Gage  grids  are  0.635  cm  square  with  ten  elements, 
each  0.005  cm  thick  and  0.05  cm  wide. 
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FIGURE  14  SCHEMATIC  DIAGRAM  OF  PLEXIGLAS-PLEXIGLAS  DOUBLE  SHOCK  EXPERIMENTS 
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FIGURE  15  SCHEMATIC  DIAGRAM  OF  WAVE  INTERACTIONS  IN  DOUBLE  SHOCK  EXPERIMENTS 
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FIGURE  16  SCHEMATIC  DIAGRAM  OF  IMPEDANCE  MATCHING  TECHNIQUE 
FOR  DETERMINATION  OF  STRESS  ASSOCIATED  WITH  SECOND 
COMPRESSIVE  PULSE 
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FIGURE  18 


OSCILLOSCOPE  RECORDS  FROM  SHOT  S7 


resistance  change  associated  with  the  second  loading.  We  base  this  on 
our  uncertainty  about  the  reloading  behavior  of  both  Plexiglas  and 
\ t terbium,  especially  around  the  yield  point  of  ytterbium.  The  magnitude 
o(  tins  uncertainty  will  be  discussed  in  Section  1.7.5.  We  expected 
approximately  8 -sec  of  one-dimensional  flow  time  in  the  gages  and  we 
believe  the  data  are  good  out  to  at  least  this  time. 

Two  other  double-shock  experiments  (S29  and  S30)  of  a different  type 
were  tried  to  see  if  there  was  a significant  amount  oT  path  dependence 
during  loading.  Figure  19  is  a schematic  diagram  of  the  target  and 
projectile  arrangement  used  to  produce  waveforms  of  the  type  shown  in 
figure  20. 

The  targets  contained  a total  of  four  gage  packages,  two  in  Type  II 
t\A  Plexiglas  and  two  in  soda-lime  glass  (impedance  = Id. 7 ± 0.1  g . mm/cm^ . 
bach  target  consisted  of  two  15-cm  diameter  0.158-cm  thick  Plexiglas  disks 
and  two  15-cm  diameter,  0.317-cm  thick  soda-lime  glass  disks.  Two  gage 
packages  were  placed  between  the  Plexiglas  disks  and  two  between  the 
glass  disks,  and  all  four  disks  were  then  cemented  together  with  epoxy- 
resin. The  leads  were  brought  out  the  sides  of  the  target. 

The  projectiles  were  fabricated  from  6061-T6  aluminum  cylinders 
10. la  cm  in  diameter  and  18  cm  long,  and  weighed  approximately  2 kg. 

The  heads  were  9.9-cm  diameter,  1.27  cm  th;  k Plexiglas  disks,  which 
were  glued  directly  into  a machined  space  in  he  front  of  the  projectile 
and  butted  directly  up  against  the  aluminum  so  that  all  stress  relief 
originated  at  the  rear  free  surface  of  the  targets.  Triggering  of  the 
power  supply  and  oscilloscopes  was  accomplished  in  the  same  way  as  in 
the  low  stress  step-release  shots.  We  expected  one-dimensional  strain 
states  in  the  target  area  occupied  by  the  gage  packages  of  about  8 .see. 


x sec ) . 


O-RING  GROOVES 


PROJECTILE:  TYPE  II  UVA  PLEXIGLAS  HEAD,  1.27  cm  THICK 
6061-T6  ALUMINUM  BODY 


PLEXIGLAS 


TARGET:  TYPE  II  UVA  PLEXIGLAS  AND  SODA-LIME  GLASS.  DIAMETER:  15  cm 

Fiberglass-epoxy  gage  package  is  approximately  2 cm  square  and 
0.025  cm  thick.  Gage  grids  Die  0.635  cm  square  with  ten  elements, 
each  0.005  cm  thick  and  0.05  cm  wide. 

MA-1797-10 

FIGURE  19  SCHEMATIC  DIAGRAM  OF  STEP-LOADING  EXPERIMENTS 
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FIGURE  20  SCHEMATIC  DIAGRAM  OF  SURFACE  INTERACTIONS  FOR  EXPECTED 
STRESS-TIME  HISTORIES  FOR  SHOTS  S29  AND  S30 


1.4.5  Phase  Transition  Experiment — A number  of  investigators 
(Hall  and  Merrill,  1963;  Hall  et  al  . 1963)  have  reported  that  the 
pressure-induced  ICC  • BCC  phase  transition  seen  in  ytterbium  is  quite 
slow  under  hydrostatic  pressure  conditions,  sometimes  requiring  many 
hours  for  completion.  This  raised  the  possibility  that  ytterbium  might 
be  usable  as  a stress  transducer  above  the  phase  transition  point  if  the 
transition  did  not  take  place  during  dynamic  loading  in  the  microsecond 
or  millisecond  time  scale.  If  the  metastable  1'CC  phase  could  be  retained 
wc  could,  theoretically,  calibrate  above  the  equilibrium  transition  point 
and  continue  to  use  ytterbium  as  a transducer  material  at  stresses  as 
high  as  we  could  go  without  causing  the  transformation.  In  addition, 
an  experiment  where  the  peak  stress  reached  was  above  the  — 39-kbar 
hydrostatic  transition  point  was  necessary  to  fix  the  upper  limit  of 
usefulness  of  ytterbium  if  the  transition  did  take  place.  Keough  (1970) 
did  one  experiment  in  which  he  double-shocked  an  ytterbium  film  to  a 
final  stress  of  approximately  14  kbar  (the  first  wave  stress  was  about 
33  kbar)  and  he  reported  seeing  a decrease  in  resistance  when  the  second 
wave  reached  the  ytterbium  and  said  that  the  decrease  continued  to  be 
gradual , not  immediate.  Wc  have  examined  Keough’s  original  oscilloscope 
records  and  found  (1)  the  resistance  decrease  attributed  to  the  phase 
change  occurs  within  300  nsec  after  the  arrival  of  ihe  second  compressive 
pulse  and  (3)  a decrease  to  about  50  per  cent  of  the  peak  resistance 
change  lcquiicd  less  than  1 ^scc.  Our  only  hesitation  in  accepting 
these  data  as  a measure  of  the  kinetics  of  the  ICC  - BCC  transition 
was  that  the  gage  showed  some  anomalous  behavior  during  the  decrease, 
and  the  shape  of  the  trace  could  have  conceivably  been  attributed  to 
gage  failure  after  the  arrival  of  the  second  compressive  pulse.  In 
addition,  Keough  assumed  that  the  resistance  change  lie  saw  at  the  shock- 
induced  phase  transition  point  corresponded  to  the  hydrostatic  phase 
change  stress,  i.e.,  - 39  kbar,  whereas  if  he  had  fit  his  data  to  a curve 
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without  that  point  and  extrapolated,  he  would  have  predicted  a shock- 

induced  phase  change  stress  of  closer  to  35  kbar,  and  the  upper  end  of 

his  curve  would  have  had  a different  shape. 

One  experiment  was  done  in  an  attempt  to  answer  some  of  the  questions 

discussed  above.  A soda-lime  glass  target  was  constructed  (impedance  - 14.7 
2 

±0.1  gm  • mm/cm  • usee)  which  contained  two  0.005-cm  thick  ytterbium  foil 

* 

grids  and  a 0.0025-cm  thick  manganin  foil  stress  gage.  It  was  impacted  with 
a glass-headed  projectile  designed  to  produce  a 1-p.sec  flat-topped  wave  of 
approximately  40-kbar  stress  intensity.  The  release  originated  at  the 
back  of  the  projectile  head,  and  triggering  was  done  with  the  usual 
combination  of  radial  pins  and  foil  switches.  The  manganin  gage  was 
included  as  a standard  for  comparison  because  it  should  produce  the 
correct  wave-form  in  a simple  material  like  glass,  and  also  should  give 
the  magnitude  of  the  stress  within  the  confidence  limits  associated  with 

manganin  gages.  The  results  of  this  experiment  are  discussed  in  Sections 
4.6.2  and  4.7.6. 

4 . 5 Static  Experiments 

Two  other  types  of  experiments  were  performed  to  ascertain  the  effects 
o-  stress  (or  strain)  on  the  electrical  resistance  of  ytterbium.  These 
were  static  (but  not  hydrostatic)  compression  tests  and  tensile  strain 
tests.  The  compression  tests  were  done  to  explore  the  possible  significance 
oi  rate  effects  in  the  relationship  of  stress  to  electrical  resistance.  If 
there  was  no  significant  difference  between  static  and  microsecond  tests, 
then  the  millisecond  response  should  be  no  different.  The  tensile  strain 
experiments  were  done  to  ascertain  the  effect  of  strain  caused  by  divergent 
flow  on  field  gage  response. 


Driver-Uarris , Inc. 
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t.J.l  at.au-  Compression  fcsts— Static  electrical  resistance  vs 
compressive  stress  tests  were  clone  on  bare  0.0025-cm  (1-mil)  foil  so 
that  results  could  be  compared  with  previous  data  (Ginsberg,  1971)  and 
on  fiberglass-encapsulated  0.005-cm  foil  to  compare  the  results  with  the 
data  I l'oin  this  program.  The  bare  specimens,  which  were  grids  identical 
to  those  used  in  the  shock  experiments,  were  glued  down  to  a a -cm  diameter 
Plexiglas  disk,  0.6  cm  thick  and  covered  with  a disk  of  tiie  same  size. 

1 h°  encapsulated  specimens  were  tested  in  the  same  manner,  but  in  one 
experiment  steel  disks  were  substituted  for  the  Plexiglas.  Pressure  was 
applied  by  a hydraulic  ram  system,  and  resistance  measurements  were  made 
with  cither  a double  Kelvin  bridge  or  Cimron  Model  6852  digital  multimeter. 

° * 2 1 (‘tis  i lo  .Strain  Tests  — Stress  waves  originating  from  the 

detonation  ot  high  explosives  or  nuclear  devices  in  or  above  the  ground 
have  a diverging  component  which  probably  causes  some  distortion  of  the 
signal  obtained  with  piezorosistant  stress  gages.  This  distortion  would 
arise  from  deformation  of  the  active  element  in  directions  contained  in 
the  plane  of  the  gage  element.  Discussions  of  this  phenomenon  are 
contained  in  Keough  (1970)  and  Kcough  et  al.  (1971),  but  these  discussions 
assume  that  the  flow  field  is  cylindrical ly  divergent  (rather  than 
spherically  divergent),  that  all  flow  in  the  gage  is  plastic,  and  that 
the  resistivity  of  the  gage  remains  constant.  Fach  of  these  assumptions 
can  have  a major  effect  on  any  analysis  of  the  real  situation.  See 
further  discussion  in  Section  1.7  and  Appendix  A. 

Our  tensile  strain  experiments  were  designed  to  be  exploratory  in 
nature;  that  is,  it  was  not  completely  obvious  which  of  many  possible 
experiments  would  be  the  most  relevant  to  the  field  gage  interpretation 
problem . fhc  configuration  finally  chosen  was  the  dog-bone  tensile 
specimen  with  attached  ytterbium  grid  and  strain  gage,  as  shown  in 
figure  21. 


DIRECTION  OF 
APPLIED  STRAIN 


FIGURE  21  DOGBONE  TENSILE  SPECIMENS 


The  dog  bones  were  mechanically  cleaned  and  then  treated  with  M-Prep 
* * 

Conditioner  and  M-Prep  Neutralizer  5 and  both  the  strain  gages  (Model 

* 

No.  KP-08-250BG-120  ) and  ytterbium  grid  were  bonded  to  the  center  of  the 

* 

gage  section  with  GA-5  or  GA-2  strain  gage  cement.  Electrical  connections 
were  made  through  copper  strips  which  were  spot-welded  to  the  ytterbium 
gage  leads.  Electrical  lead  wire  was  then  soldered  to  the  copper.  Strain 
was  applied  to  the  smaller  dog  bones  with  a screw-driven  tensile  tester, 
and  to  the  larger  ones  with  an  hydraulically  driven  apparatus. 

The  120-J  strain  gages  were  monitored  by  a Wheatstone  bridge  circuit 
driving  the  X-axis  oi'  a high-impedance  X-Y  recorder.  One  leg  of  the 
bridge  was  the  active  gage  and  the  opposite  leg  was  an  identical  gage 
mounted  on  a similar  aluminum  test  strip  to  provide  gage  temperature 
compensation . 

The  1.5  u ytterbium  gages  could  not  be  measured  with  a standard 
Wheatstone  bridge  because  of  their  low  resistance,  so  a bridge  circuit 
with  a high  gain  amplifier  was  used  to  provide  an  adequate  signal  to 
drive  the  Y-axis  of  the  X-Y  recorder.  This  arrangement  made  it  possible 
to  plot  the  change  in  resistance  of  the  ytterbium  gages  as  a function  of 
the  specimen  strain,  which  was  the  information  we  wished  to  obtain  from 
these  experiments. 

■1.3.3  Direct  Divergent  flow  Simulntion--Simulation  of  spherical 

divergent  flow  was  attempted  by  cementing  ytterbium  grids  to  the  inflatable 
t 

sphere  shown  in  figure  22  and  the  change  in  electrical  resistance  monitored 
as  a function  of  local  strain  (measured  optically)  as  the  sphere  was 
inflated.  Although  changes  were  seen  in  resistance  as  the  sphere  size  was 
changed,  we  were  unable  to  show  conclusively  that  the  strain  measured  on 
the  sphere  surface  adjacent  to  the  ytterbium  grid  was  equal  to  or  had 


* Micro-Measurements,  Inc.,  Romulus,  Michigan. 

! Sun  Products  Corp.,  Barberton,  Ohio  11203,  No.  10135. 
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some  simple  relationship  with  the  strain  in  the  grid,  and  we  were  unable 
to  measure  the  actual  grid  strain  with  sufficient  accuracy.  This  experi- 
ment was  abandoned , but  the  same  type  of  experiment  using  a stiffer  sphere 
material  might  give  useful  results. 


•l.G  Results 

This  section  presents  the  results  of  each  of  the  types  of  experiments 

performed  to  study  f he  electrical  resistance  of  ytterbium  as  a function 

ol  stress,  strain,  or  temperature.  These  results  are  discussed  in  detail 

in  the  Interpretation  and  Discussion  Section  (Section  d.7)  that  follows. 

>.6.1  Dynamic  Calibration  Experiments — The  results  of  our  dynamic 

calibration  shots  are  shown  in  Tallies  i and  5 and  figures  23  through  dl. 

Talile  I contains  all  of  the  data  points  obtained  during  the  course  of 

this  program  and  some  that  were  obtained  between  the  completion  of  our 

previous  el  fort  (Ginsberg,  1971)  and  the  beginning  of  the  present  program 

from  shots  that  were  performed  to  solve  specific  calibration  problems.  The 

latter  are  designated  as  the  IRD  series  shots.  All  data  that  were 

rejected  are  identified,  and  the  reasons  for  rejection  are  discussed 

later  in  the  report.  The  loading  data  are  the  set  of  all  points  reached 

by  single  shock  loading  to  a steady-state  end  state,  and  consist  of  the 

set  of  AR/R^  and  q values  contained  on  the  first  line  of  the  data 

identified  with  each  shot.  The  unloading  data  are,  for  each  shot,  the 

set  ol  points  (starting  with  the  peak  stress)  contained  in  the  columns 

labeled  Gage  1,  AR/R  and  Gage  2 AR/R  . Brief  information  of  particular 

o o 

interest  about  each  shot  is  included  in  the  comments  where  appropriate. 

fable  u contains  the  function  (jt  ~ i (AR/R  ) w'e  have  chosen  to  represent 

* 1 1 ° 

each  ol  the  unloading  paths  and  the  standard  deviation  in  stress  assoc- 
iated witli  each  fitted  equation. 


* 


Our 
4.7; 
s uni 


criteria  for  the  choice  of  the  function  are  discussed  in  Section 
the  choice  is  not  always  the  best  fit  in  terms  ol  minimizing  the 
of  the  squares  of  residuals  about  the  mean. 
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SUMMARY  OF  DYNAMIC  CALIBRATION  EXPERIMENTAL  DATA 


•o  re 
re  X 


O 

u 

a 


re 

x 


*o 
re  o 
W5  C 
O *H 
pH  C3 


Ifi 


c 

o 

E 

E 

0 

O 


CM 

0 

X 

pH 

[p 

X 

CM 

CM 

fp 

CM 

O 

a 

X 

co 

3* 

CM 

rH 

CD 

CD 

o 

b£ 

re 

\ 

*T 

X 

X 

rH 

0) 

CM 

01 

tp 

O 

<j 

pH 

c 

o 

d 

3* 

CM 

o 

o 

o a 

ti  V 

re  a 

O <J 


o 

a 


u 

o 

.3 

0 E 


X * 


a 

o 

o 

£ 


O 5h  0 

■m  re  3 

O £ H 

* a re 

> 


pH  a 

• O 
pH  0 
.C 


T3  H re 


W 

re 

o 

o 

X 

w 


O 

a 


0 

13  O 
re  (D  rH 
0 3 fc£ 

-H  pH  G 

<D  3 H 
a > re 


•o 

re 

o 


0 

a 


x cm 
00  T 
x cm 
o o 


CM  O (D  O) 
T O h-  o 
o*  x h h 
o o o o 


o o o o o o o 


CM 

i 

X 

01 

r-, 

pH 

Si 

rH 

X 

X 

i 

tp 

rH 

CM 

o 

31 

Jp 

CM 

X 

oi 

i 

pH 

o 

CP 

X 

X 

X 

X 

o 

o 

i 

i 

o 

rH 

o 

o 

o 

o 

O 

o 

o 

i 

o 

o 

o 

o 

o 

o 

o 

o 

O)  lO  00  CM 
CM  3*  X h* 
T CM  H 
O O O O 


o o o o 


CO 

pH  r» 

00  co 


CM  O 01  3*  X 

co  «5  a)  o)  m 

CO  O (p  IT3  pH 


ao  co 

CO  CM 

o o 


X CD  00  Ol 
O in  Cl  pH 
X CO  pH  H 

o o o o 


tp  N 
CM  00 
3*  pH 

o o 


P ° o O t n o o o o o o o o’  o’  o o’  o’  o’ 


Ip  pH 

CM  3 
X CO 

o o 


* 

CM  CO  CM 
CD  H in 
X co 
o o 


o 

o 


Is-  . , 

00 


00  CD 

O {p  __  _ 
Qo  ire  co  cm 
o o o o 


oooo  ooo  oooo 


re 

re 

CM 

3*  3* 

(P 

X 

o 

IP 

CM 

* 

DM 

0 

u 

u 

re 

X 

rp  cm 

01 

-** 

Ol 

o 

X 

H 

X 

X 

7 

-1 

0 

01 

fp 

Ol’ 

5 

0 

rH 

d O 

X 

a 

pH 

CM 

pH 

pH  O 
pH  fp 


*r  x 

3*  CM 


I/O  pH 
N pH 


X O X X 
Ol  W h-  *T 


*p  CO 
CM  X 


CD  pH 

3 o m re 


pH  O O O rH 


O pH  pH  O O pH  O O pHOOO 


X 

o 

I o 

w re  x 

re  re  re 


re 


re 


tL  t£  O 


re  o 

re  rH 
rH  0 

t£  U 


X 

V 

\ 0 
re  re  x 
re  re  re 


b£  b£  <D 


O re  a 

3 3*-- 

J O’  W 


J O'  X 


X X JD 
0 0 3 

pH  pH  0 

a a q 


w x a 
re  o o 


x x .Q 

0 0 3 


o a x 


pH  pH  0 

a a Q 


re 
re  x 
re  o 


o 

re 

re 

re  o 

re  pH 
pH  o 
bL  U 


re  re 
re  re 


x x xi 

0 0 3 


O a x 


pH  rH  O 

a a a 


o a x 


CO 

X 


X 

X 


CD 

X 


tp 

X 


X 

X 


63 


iaiilMlitfliiiaiil 


Mli 


NONLINEAR  UNLOADING  PATHS 
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FIGURE  23  DYNAMIC  LOADING  DATA 

Dashed  line  is  fit  to  data  from  this  work. 
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FIGURE  24  COMPARISON  OF  PREDICTED  DYNAMIC  BEHAVIOR 
WITH  EXPERIMENTAL  DATA 
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FIGURE  27  LOADING  AND  RELEASE  PATHS,  SHOT  SI 
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FIGURE  28  LOADING  AND  RELEASE  PATHS,  SHOT  S18 
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FIGURE  29  LOADING  AND  RELEASE  PATHS,  SHOT  S19 
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FIGURE  31  LOADING  AND  RELEASE  PATHS,  SHOT  S2 
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FIGURE  32  LOADING  AND  RELEASE  PATHS.  SHOT  S17 
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FIGUHfc  J5  COMPARISON  OF  HYDROSTATIC  AND  DYNAMIC  BEHAVIOR  BELOW  10  kbar 
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FIGURE  37  LOADIN  HATH  AND  SOME  REPRESENTATIVE  MEASURED  RELEASE  PATHS 
(10  ' u > 3 kbar)  0.005-cm  (2-mil)  FOIL 

Other  release  paths  have  been  measured  in  this  stress  range,  but  have  be  ;n  omitted 
for  clarity, 
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FIGURE  38  SLOPE  OF  LEAST-SQUARES  FIT  TO  LINEAR  UNLOADING  DATA 
AS  A FUNCTION  OF  PEAK  STRESS 
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FIGURE  39  FAMILY  OF  UNLOADING  PATHS  FOR  PEAK  STRESS  OF 
10  kbar  OR  LESS 
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FIGURE  40  ZERO-STRESS  RESIDUAL  RESISTANCE  (AR  /R  ) VERSUS  PEAK  STRESS 

no 
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1'iRUre  23  contains  the  curve  chosen  to  represent  the  full  set  of 
data  up  to  about  33  kbar  and  contains  our  data  points  as  well  as  those 
ol  some  other  investigators  (Keough,  1970;  Brown,  1970).  Figure  24 
compares  our  dynamic  data  with  the  dynamic  piezoresistance  function 
LUley  and  Stephens  (1971)  calculated  from  their  hydrostatic  data  by 
making  corrections  for  temperature  and  geometry  effects  and  Figure  25 
shows  our  data  and  Lilley  and  Stephens'  hydrostatic  data,  as  does 
Figure  35,  but  at  lower  stresses.  Figures  26  through  32  are  the  unloading 
data  from  peak  stress  levels  above  10  kbar  along  with  the  fitted  curves. 
Figure  33  contains  the  loading  data  for  peak  stress  values  up  to  10  kbar, 
and  Figure  31  has  a comparison  between  these  data  and  the  data  we 
obtained  previously  on  0.0025  cm  (1-mil)  foil  (Ginsberg,  1971).  Figures 
db,  3 1 , 38,  and  39  show  unloading  data  for  a < 10  kbar.  Figure  -10 
contains  t lie  values  of  the  zero-stress  residual  resistance  (AH  /R  ) as 
a function  of  peak  stress,  and  Figure  41  presents  these  vnlurs^a 
different  manner,  showing  the  zero-stress  residual  resistance  divided 
by  the  total  resistance  change  at  peak  stress.  Figures  42,  43,  and  44 

contain  the  results  of  our  comparison  experiments  between  annealed  and 
cold-worked  foil. 

4.6.2  ~iasc  Transition  Experiment— Figure  45  shows  the  records 
obtained  during  Shot  S20  which  was  fired  to  explore  the  kinetics  of 
the  high  pressure  phase  transition  seen  during  static  experiments  and 
to  ascertain  the  stress  level  at  which  it  .tarts.  Figures  45a  an.  45b 
are  the  ytterbium  gage  records,  while  45c  is  the  record  from  the  manganin 
Sage  in  the  same  plane  ns  the  two  ytterbium  gages.  Figure  46  is  a combined 

Plot  of  the  data  and  shows  some  of  Ihe  more  interesting  features  of  the 
experimental  records. 

1.6.3  Lp.n*ile  Strain  Experiments— Figures  47  through  49  contain 
the  results  of  the  experiments  we  performed  to  obtain  data  applicable 

to  the  solution  of  the  divergent  flow  problems  encountered  in  field  gage  use, 


87 


LOADING  CALIBRATION  CURVE 
0 005-cm  (2-mill  FOIL  AS-ROLLED 


NOTES:  Solid  lines  with  arrows  are  unloading  paths. 

Error  bar  represents  - one  standard  deviation 


MA-1 797-48 


FIGURE  44  COMPARISON  OF  RESPONSE  OF  ANNEALED  AND  AS-ROLLED 
(COLD-WORKED)  YTTERBIUM  FOIL 
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lal  YTTERBIUM  GAGE  NO  1. 
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(b)  YTTERBIUM  GAGE  NO  2 

Vertical  sensitivity  = 0 35  V/cm 
Horizontal  sweep  = 0.5  us/cm 


Id  MANGANIN  GAGE 
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FIGURE  45  OSCILLOGRAMS  FORM  SHOT  S20 

Peak  stress  ' 40  kbar  All  gages  in  same  plane  in  Target 
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CIGURE  46  PLOT  OF  DATA  FROM  SHOT  S20  SHOWING  PHASE  TRANSITION 
IN  YTTERBIUM 
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FIGURE  48  FRACTIONAL  CHANGE  IN  RESISTANCE  VERSUS  TENSILE  STRAIN. 
PERPENDICULAR  AND  PARALLEL  GRID  ORIENTATION 


GAGE  SECTION 


FIGURE  49  FRACTIONAL  CHANGE  IN  RESISTANCE  VERSUS  TENSILE  STRAIN, 
ANNEALED  AND  AS-ROLLED  YTTERBIUM 


Figure  17  shows  the  results  obtained  using  three  different  gage  section 
widths  and  also  contains  some  unloading  paths.  ligure  49  compares  the 
behavior  of  annealed  and  cold-worked  (as -rolled)  foil,  and  Figure  48 
compares  the  response  of  a parallel-mounted  grid  to  a grid  bonded  with 
its  elements  perpendicular  to  the  direction  of  tensile  strain. 

■1.6.4  Static  Compression  Piezoresistance  Tests — Figure  50  contains 
the  results  of  the  static  compression  tests  done  to  explore  the  possibility 
of  the  presence  o*  rate  effects  in  the  pic/oresistant  response  of  ytterbium. 
Figure  51  contains  some  results  obtained  by  Smith  et  al . (1972)  for 
comparison  with  our  data. 

1.7  Discussion  and  Interpretation 

In  this  section  the  results  obtained  in  this  program  are  analyzed 
and  discussed.  First  the  aspects  of  the  analysis  that  were  common  to 
most  of  the  experiments  are  considered,  such  as  curve-fitting  techniques 
and  estimation  of  experimental  errors.  Each  of  the  specific  types 
of  experiments  performed  arc  treated,  and  finally  some  general  comments 
are  given  about  thi  piezoresistance  of  ytterbium. 

1.7.1  Uncertainties  in  Experiments  and  Annlysis--Two  measured 
quantities  are  important  to  a calibration  program.  These  are  the  measure- 
ment of  the  ratio  of  change  in  resistance  to  the  original  resistance  of 

the  ytterbium  gage,  or  ilt/H  , and  the  stress  levels  at  which  we  wish  to 

o 

know  AR/R  . The  uncertainties  in  the  measured  value  of  AR/R  can  arise 
o o 

from  (1)  possible  systematic  or  random  electronic  equipment  malfunctions, 

(2)  oscilloscope  calibration  errors,  (3)  measurement  of  circuit  parameters 

* 

in  the  AV/V  to  AR/R  correction,  (4)  conversion  of  the  oscillograms  to 
o o 

digital  form,  and  (5)  uncertainty  in  the  estimation  of  the  best  value 


of  AV/V  for  a given  value  of  stress  from  the  plotted  peak  stress  and 
o 

release  step  DC  levels. 


* See  Appendix  B. 
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kbar 


FIGURE  50  STATIC  PIEZORESISTANCE  DATA  COMPARED  WITH  DYNAMIC  DATA 
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STRESS  — kbars 


GA  8567-75 


FIGURE  01  STATIC  CALIBRATION  CYCLES  FOR 
AN  YTTERBIUM  STRESS  GAGE, 
from  Smith,  et  al.  (1972) 
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tte  must  assume  that  we  can  recognize  the  occurrence  of  any  syste- 
matic or  random  electronic  malfunctions  that  could  be  large  enough  to 
cause  major  errors.  Our  systems  and  experiments  had  sufficient 

redundancy  and  cross-checking  capability  to  eliminate  this  possibility. 
Oscilloscope  calibration  errors  can  occur  if  there  is  some  important 

change  in  the  recording  system  between  the  time  of  the  experiment  and 

the  time  of  calibration,  and  if  there  are  nonlinearities  in  the  system 

that  are  not  taken  into  account.  Wo  tried  to  do  our  calibrations 

immediately  after  each  shot,  and  also  checked  each  set  of  calibrations 

for  nonlinearity  from  the  top  to  the  bottom  of  the  trace.  Occasional 

dii terences  wore  noticed  of  up  to±2  , between  readings  from  the  same  gage 

on  different  oscilloscopes  and  some  changes  in  calibration  from  shot  to 

shot  on  the  order  of  ±1  . Measurement  of  circuit  parameters  used  in  the 

conversion  of  l\'/\  values  to  the  actual  Mi  R values  was  to  better  than 

o o 

±1  . Conversion  of  the  raw  data  to  digital  form  adds  practically  no  error. 

" e c'st imate  the  value  of  IV  lor  each  level  by  making  a number  of 

readings  at  different  times  and  averaging  the  result.  The  average 

deviation  around  the  mean  LX  V for  each  level  was  usually  less  thin  ±1T. 

there  is  no  standard  way  to  add  up  these  uncertainties,  so  we  will 

claim  that  the  accuracy  in  any  single  measured  value  of  £R/R  at  any 

o 

given  experimental  stress  level  was  between  ±2r;  and  ±r>r’r,  with  the  higher 
value  applicable  to  the  more  uncertain  cases,  sucli  as  the  last  step  in 
a step-release  experiment,  or  a zero-stress  residual  resistance  value. 

1 he  identifiable  errors  in  stress  in  the  dynamic  compression  experi- 
ments arise  from  uncertainties  in  the  measured  values  of  projectile 
velocity  and  subsequent  conversion  to  particle  velocity  and  uncertainties 
in  the  equations  of  state  used  for  the  target  and  flyer  materials  through 
which  the  particle  velocity  is  related  to  stress.  We  can  measure  the 
projectile  velocity  with  an  accuracy  of  about  ±2^.  Barker  and  Ilollenbach 
(1970)  claim  an  accuracy  of  ±2%  in  their  plexiglas  and  quartz  equation 
of -state  data  and  our  glass  data  are  accurate  to  this  level  at  the  low 
stresses  encountered  here.  We  expect,  therefore,  that  an  accuracy  of 
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±)'c  in  stress  is  reasonable  for  each  of  the  measured  stress  levels  in 
the  calibration  experiments. 


*•7.-  ( urve  1 1 1 tl ng  Each  of  the  dynamic  calibration  experiments 

was  designed  to  give  at  least  two  points  on  the  loading  curve,  while  the 

majority  of  the  shots  were  designed  to  produce  unloading  data  as  well. 

A typical  low-stress  step-release  experiment  (glass  flyer,  Plexiglas  target) 

gave  two  loading  points  and  eight  unloading  points.  Our  ultimate  choice 

of  treatment  of  the  data  was  based  on  a combination  of  statistical  criteria, 

convenience  of  computation,  and  physical  reasonableness.  Note  that  the 

loading  curves  are  given  in  the  form  n = f(AK/R  ) rather  than  the  reverse 

o ’ 

whereas  AH  lt^  is  actually  a function  of  stress.  This  was  done  to  simplify 

the  use  of  the  analytical  forms  in  reading  stress  gage  outputs.  If  we 

were  to  try  to  look  for  a physical  model  to  explain  the  piezoresistance 

of  ytterbium  we  would  refit  the  curves  and  put  them  in  the  form  AR/R  = f (q-) 

o 

The  unloading  curves  were  fitted  to  straight  lines  for  peak  stresses 
up  to  10  kbar  by  inspection  without  using  any  "best  fit'  criteria,  i.e., 
we  did  not  try  any  other  functions.  It  is,  however,  obvious  that  the 
measured  paths  are  linear  or  are  very  close  to  being  linear.  Fitting  the 
unloading  paths  originating  at  peak  stress  values  above  10  kbar  was  more  of 
a problem;  this  will  be  discussed  in  the  next  section. 

In  the  discussion  that  follows,  we  use  the  terms  experimental  value  of 
stress  to  signify  the  stress  calculated  from  the  known  properties  of  the 
target  and  flyer,  and  the  measured  projectile  velocity.  First,  the  £.R/R 

o 

values  are  converted  to  stress  using  one  of  the  calibration  equations 
(equation  1 or  2) . Then,  the  stress  derived  in  this  way  is  compared  to  the 
experimental  stress.  The  difference  is  an  indication  of  the  accuracy  of 
each  data  point . 

1 •7*3  High-Stress  Step-Release  Experiments  — Shots  SI  and  S2  were 
our  first  lucalox  flyer-quartz  target  step-release  experiments.  We 
reached  a peak  stress  of  13.82  kbar  during  SI  and  29.97  kbar  during  S2. 

The  data  from  SI  were  not  of  the  expected  quality  because  a substantial 
amount  of  shock  tilt,  coupled  with  the  effects  of  the  relatively  thick 
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^ubout  0.020  err.)  gage  package,  marie  the  steps  less  sharp  than  we  would 


have  liked.  The  peak  'UR  values  lor  the  two  gages  were  1.192  and 

o 

1.185,  respectively,  the  rounded  values  ol'  which  coi  respond  to  a peak 

stress  of  13.38  kbar,  a deviation  of  -0.11  kbar  from  the  experimental 

stress  of  13.82  kbar.  The  release  data  plotted  in  Figure  27  also  show 

good  reproducibility.  he  had  corrected  the  tilt  problem  by  the  time  we 

fired  S2,  so  the  data  were  more  easily  read  from  the  oscillograms.  The 

rounded  peak  values  of  HR  for  the  two  gages  were  1.33  and  1.91, 

o 

corresponding  to  deviations  of  -1.11  and  1.72  kbar,  respectively,  or 
l.tTo  and  G.l'l  ol  the  experimental  value  of  29.97  kbar.  The  release 
data  obtained  during  S2  and  shown  in  Figure  31  were  somewhat 
difficult  to  represent  by  a polynomial  that  fit  well  and  did  not  have 
too  many  inflection  points,  he  decided  that  the  unloading  curve  would 
be  lit  reasonably  well  by  the  loading  curve  but  with  a zero-stress  offset, 
so  we  added  a number  of  the  calculated  loading  points  to  the  measured 
population  of  unloading  data,  and  weighted  the  zero  stress  value  ol  /R/R 

o ’ 

Ihe  fitted  line  crosses  the  loading  path  twice,  but  it  is  a reasonable 

representation  of  the  unloading  behavior  from  this  stress  range. 

The  rest  of  the  series  of  Lucalox-quartz  step  release  shots  were 

•’ll  and  812,  and  SI/  through  S22  (excluding  S20).  We  obtained  acceptable 

loading  and  release  data  from  both  gages  in  all  shots  except  S17,  in 

which  the  peak  stress  was  33.10  kbar.  One  of  the  gages  in  S17  behaved 

in  the  expected  manner,  but  the  other  showed  a rise  to  a relatively  low- 

change  in  resistance,  then  a sudden  drop  to  a lower  plateau,  where  it 

stayed  for  approximately  one  „sec  until  the  arrival  of  the  first  release 

step.  We  rejected  this  value  because  we  believed  that  the  ytterbium 

making  up  this  gage  had  begun  to  transform  into  the  high  pressure  (HCC) 

phase  and  was  showing  the  behavior  noted  when  the  phase  change  in  S20 

was  studied.  It  is  interesting  to  note  that  the  other  gage  din  ^ot  seem 

to  begin  to  transform.  We  inferred  that  the  phase  change  probably 

nucleates  at  local  stress  concentrations,  such  as  around  inclusions,  and 
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quickly  grows  from  those 
uncertainty  in  predicting 
uill  begin,  and  leads  to 
stress  at  whicii  commercia 


nuclei.  It  is  this  phenomenon  which  leads  to  our 
the  macroscopic  stress  at  which  the  transformation 
our  recommendation  that  30  kbar  be  the  maximum 
lly  available  ytterbium  should  be  used  without 


a backup  system. 

The  best  accuracy  was  obtained  in  S21 , where  the  derived  peak  stress 
values  of  the  two  gages  deviated  by  2.8%  and  2.2";  from  the  experimental  value, 
and  in  S17  where  the  deviation  was  2.1%  for  the  gage  that  functioned.  The 
poorest  accuracy  was  exhibited  by  the  data  from  S22;  here  the  derived  value 
of  stress  from  LH  Hq  of  the  first  gage  deviated  by  11.4%  from  the  experimental 
value;  we  did  not,  however,  have  any  grounds  for  rejection  of  these  data. 

Shots  SI 1 and  S12  were  done  as  a pair  to  chock  (as  well  as  could  be 
done  with  one  pair  of  shots)  the  effect  of  the  fiberglass-epoxy  package 
on  the  response.  The  target  for  Sll  contained  two  bare  ytterbium  grids; 
that  for  S12  two  of  our  regular  liberglass-epoxy  encapsulated  grids.  We 
tned  for  25  kbar  in  Sll  and  obtained  2 1.03,  which  was  then  the  stress 
aimed  for  in  S12,  where  we  obtained  23.81  kbar.  Given  our  ±3%  estimate 
of  accuracy  in  stress,  the  two  can  bo  considered  to  be  the  same  for  purposes 
of  comparison.  The  values  from  the  two  encapsulated  gages  were  both 

lower  than  either  of  the  readings  from  the  bare  gages,  but  since  the 
highest  reading  showed  a considerable  deviation  from  the  mean  and  the 

other  three  formed  a cluster,  no  difference  in  behavior  could  be  inferred 
from  this  pair  of  shots. 

There  is  one  other  interesting  aspect  of  the  Lucalox-quartz 
experiments.  On  examining  the  records  from  SI,  we  noticed  that  the 
recording  time  (over  8 use,.'  was  significantly  longer  than  expected, 
and  we  also  noticed  this  in  subsequent  shots  of  this  type.  We  had 
expected  to  he  able  to  obtain  a value  of  AK/Hq  corresponding  to  the 
peak  stress  and  three  release  calibration  steps  while  the  target  and 
enough  of  the  projectile  head  were  in  one-dimensional  strain  states. 

This  was  done  in  all  experiments  except  SI,  where  dispersion  arising 
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I rom  tilt  obscured  the  third  step.  However,  a zero-stress  residu.il 

resistance  ( MI  /R  ) measurement  was  not  expected  on  these  shots 
11  o 

because  the  targets  were  not  expected  to  last  long  enough.  Given  the 

long  recording  time,  we  did  measure  a value  of  £11  R at  a stress  level 

o 

very  close  to  zero  and  included  it  with  the  other  ill  R data.  Two  things 

II  o 

must  be  remembered  here:  the  stress  was  not  actually  zero  (but  is  always 
under  1 kbar) , and  the  target  and  gages  were  not  in  one-dimensional  strain. 
However,  since  the  gages  showed  no  obviously  suspicious  behavior,  we 
decided  to  include  these  data  with  the  rest  of  the  zero-stress  residual 
resistance  data. 

Tsually  the  unloading  paths  were  easily  fit  by  polynomials  in  Ill'll  . 

o 

In  some  cases,  however,  fit  was  sacrificed  for  physical  reasonableness. 

* 

High  order  polynomials  often  fit  the  data  well,  but,  given  the 
relatively  small  number  of  points,  tend  to  fit  very  poorly  between 
the  experimental  data  and  not  to  be  physically  reasonable  because 
of  the  large  number  ol  inflection  points.  The  only  release  path  that  had 
to  be  force-fit  in  any  way  was  that  from  Shot  S2,  as  explained  previously. 

The  standard  deviations  in  stress  associated  with  the  fitted  release 
paths  vary  from  0.2(15  kbar  in  S17  to  0.8G5  kbar  in  S12.  One  way  of  treating 
this  information  would  he  to  convert  this  value  to  an  indication  of  the 
percentage  uncertainty  at  the  midpoint  stress — 1G.7  kbar  in  S17  and  15  kbar 
in  S12 . The  percentage  uncertainties  in  stress  at  these  midpoint  stress 
levels  are  then  1.3 .Vi  and  5.7G  , respectively. 

In  summary,  we  were  successful  in  using  the  I.ucalox-quartz  step 
release  experiments  to  produce  a large  amount  of  loading  and  release 
data  in  the  range  of  8 kbar  and  to  over  33  kbar  despite  the  resolution 
difficulties  caused  by  the  presence  of  the  thick  gage  package  and  the 
transient  disturbances  to  the  flow  that  were  present  during  the  experiments. 


* 

Using  a least-square  criterion. 
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'•7<*  Low-Stress  Step-Release  Experiments — Another  step-release 
configuration  (glass  i lyer-Plexiglas  target)  was  used  to  obtain  all  of 
our  release  data  and  a large  part  of  our  loading  data  at  stiess  levels 
of  6 kbar  and  below.  The  first  series  consisted  of  Shots  S4  , S6,  S8 , S10, 
S13,  SM,  and  S15.  Accuracies  ranged  from  excellent  in  Sid,  where  the  two 
gage  readings  deviated  by  0.08  kbar  (0.17%)  and  0.02  kbar  (0.04%)  from 
the  experimental  value  of  4.72  kbar  and  S6  and  S4,  where  the 
accuracy  was  of  the  same  order,  to  poor  in  one  shot  - S10.  In  this 
case  the  percentage  deviations  from  the  experimental  stress  were  16.5%  and 
13.1  lor  the  two  gages.  We  believe  (and  will  discuss  this  in  more  detail 
later)  that  discrepancies  of  this  type  arise  from  differences  in  the 
metallurgical  history  of  the  foil  used  to  make  the  gages.  However, 
because  we  had  no  practical  and  reliable  way  of  testing  each  gage  for 
some  indication  of  this  before  each  shot,  we  can  only  suggest  this  as  a 
reason  lor  the  occasional  lack  of  reproducibility. 

The  series  was  concluded  with  Shots  S21  and  S25,  which  were  done  to 
check  for  possible  time  effects  by  changing  the  flyer  thickness  and 
thereby  the  step  length,  and  S26  through  S28,  which  were  comparisons 
between  annealed  and  as-rolled  ytterbium  foil.  The  only  data  that  were 
i ejected  were  those  from  Gage  2 in  S13,  where  the  derived  stress  value 
was  found  to  differ  from  the  expected  value  by  much  more  than  two  standard 
deviations,  a criterion  for  rejection. 

Our  treatment  of  the  unloading  data  derived  from  these  experiments 
will  be  discussed  first,  hach  shot  (except  S13  and  the  annealed  versus 
cold-worked  shots)  produced  eight  points  that  comprised  the  data  for  one 
release  path  originating  at  the  peak  stress  reached  in  the  experiment. 
These  points  were  plotted,  and  straight  lines  were  fit  to  all  of  the 
paths  up  to  and  including  the  Lucalox-quartz  Shot  S22  (8.22  kbar).  Some 
of  the  results  are  shown  in  Figures  36  and  37.  We  now  had  a set  of  14 
straight  lines  (including  Shot  XI  — see  Table  4)  of  the  form 
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where  a Is  the  stress  along  the  release  path,  A and  li  are  constants  lor 

each  individual  path  hut  are  1 mictions  r>!  the  parameter  a the  peak 

!>’ 

stress  reached  in  the  experiment,  and  (I  1 > i-  the  fractional  resistance 


change  along  the  unloading  path  « Vti  t )i«  *n  }>  1 
(Figure  38)  and  found  the  best  titling  1 
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To  find  the  best  values  <il  \,  *.•  noted  that 


ilH 

II 


A 

11 


(11) 


toi  anj  Hi  veil  peak  stress  q , where  I!  II  is  the  zero-stress  residual 

P II  > 

resistance.  Our  measured  J.K  li  values  wore  lit  to  a number  of  functions 

H o 

and  an  acceptable  fit  was  found  (see  Figure  |o)  to  be 
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Unfortunately,  the  problem  is  now  over-defined.  We  have  a value  for 
the  slope  and  both  end  points  of  every  line  as  a function  of  peak  stress, 
and  if  we  use  any  pair,  the  value  of  third  quantity  does  not  fit.  The  only 
reasonable  answer  we  can  Rive  at  this  time,  pending:  further  analytical 
work  on  the  data,  is  to  say  that  it  is  most  important  that  the  unloading 
paths  go  through  the  loading  curve  at  the  proper  peak  stress,  and  the 
zero-stress  end  point  or  slope  values  must  conform  to  this. 

We  have  constructed  a set  of  unloading  paths  by  drawing  straight  lines 
through  the  zero-stress  ond  points  and  the  corresponding  peak  stress  values. 
This  family  of  lines  is  shown  in  Figure  39.  Hearst  (1971)  has  carefully 
examined  the  data  contained  in  this  report  and  constructed  a set  of  algorithms 
for  computerized  reduction  of  gage  data  in  both  loading  and  unloading.  He 

reports  that  they  have  been  used  and  are  satisfactory. 
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The  release  data  presented  in  our  previous  report  (Ginsberg,  1971)  which 
were  tentative,  should  not  be  used.  An  error  arose  from  our  inability  to 

describe  the  equation  of  state  of  Plexiglas  in  a zone  where  there  was 

* 

interaction  between  compression  and  release  wavcu,  thereby  causing 
both  time  and  stress  level  shifts  in  our  code  output.  The  release  paths 
presented  here  were  derived  from  good  experimental  data  and  can  be 
used  to  obtain  stress-time  unloading  data  within  the  stated  limits  of 
accuracy. 

The  available  recording  time  in  this  series  of  shots  was  much  longer 
than  the  4 p,sec  of  one-dimensional  strain  time  in  the  projectile  head. 

The  data  obtained  from  steps  as  long  as  7 (j,sec  from  impact  were  consistent 
with  the  data  taken  at  earlier  times.  It  should  be  pointed  out,  however, 
that  the  gages  themselves  had  approximately  6 ysec  of  time  in  one-dimensional 
strain.  Shots  S24  and  S25  were  done  at  about  the  same  stress,  but  with 
two  different  flyer  plate  thicknesses,  0.317  (S24)  and  0.635  cm  (S25). 

After  about  4.5  Usee,  the  projectile  heads  were  both  out  of  the  one- 
dimensional strain  time  region,  but  S24  had  undergone  four  release  steps 
whereas  S25  had  undergone  only  two.  The  rest  of  the  steps  in  S25  took 
place  under  mixed  strain  conditions  in  the  projectile  head.  However, 
no  apparent  difference  was  seen  in  the  behavior  of  the  gages  in  the  two 
shots  under  these  two  different  conditions. 

In  our  previous  report  (Ginsberg,  1971)  we  suggested  that  the 
metallurgical  history  of  the  foil  used  in  the  experiments  was  likely  to 
be  important  in  that  variations  in  history  would  probably  cause  variations 
in  electrical  response  to  stress,  thereby  compromising  both  the  accuracy 
of  the  gages  and  the  precision  of  the  calibration.  Such  differences  have 
been  shown  to  be  a problem  with  manganin  (Keough  and  Wong,  1970;  Rosenberg 
and  Ginsberg,  1972)  and  are  not  really  surprising  nor  unexpected.  Shots 
S26 , S27 , and  S28  were  designed  to  examine  the  magnitude  of  the  difference 

■ ■ ■■  — ■ AD 

* Our  analysis  was  based  on  going  from  — - (t)  to  a(t)  by  computer  code 

R 

o 

predictions,  using  available  equation  of  state  data. 
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in  dynamic  piezoresistant  response  between  as-rolled  foil  (which  is  the 

usual  condition)  and  foil  which  had  been  annealed  at  approximately  175°C * 

-5  t 

lor  30  minutes  at  a pressure  of  10  mm  Hg.  The  results  showed  that  there 
is  a very  large  difference  in  the  behavio'  of  these  foils.  Plots 
of  two  of  the  oscillograms  are  contained  in  Figures  42  and  43,  and  an 
approximate  loading  curve  for  annealed  foil,  along  with  measured  release 
paths,  is  shown  in  Figur  44.  Vie  had  seen  previously  that  the  response 
ol  the  as  rolled  AREF  fc.il  is  close  to  the  same  as  that  of  the 
UC  foil,  so  the  difference  can  be  taken  to  be  between  annealed  and 
cold-worked  foil.  The  magnitude  of  the  difference  in  response  suggests 
that  the  electrical  resistivity  of  ytterbium  is  very  sensitive  to  the 
concentration  of  strain-induced  lattice  defects. 

i he  data  from  Shot  S28  also  (Figure  43)  are  of  interest  here.  Note 
,l*at  although  the  tops  of  both  outputs  are  distorted  by  a spike,  that  of 
the  rolled  foil  gage  is  relatively  flat  while  the  other  is  climbing  at  a 
rapid  rate.  We  would  like  to  suggest  the  possibility  that  defect  multi- 
plication is  taking  place  during  the  two-p,sec  intervals  when  the  gages 
are  at  the  peak  stress  of  1.40  kbar,  and  this  effect  is  causing  the 

resistance  change.  We  believe,  therefore,  that  the  value  of  AR/R  for 

o 

the  annealed  gage  is  not  an  equilibrium  value  and,  if  it  were,  the 

dii ierence  between  it  and  that  of  rolled  foil  gage  would  be  much  greater. 

The  ratio  of  AR/R  of  the  annealed  to  the  rolled  foil  in  the  three  shots  is: 
o 

at  5.36  kbar,  2.45;  at  3.21  kbar,  3.58;  and  at  1.40  kbar,  1.37.  We  must 

conclude  that  variations  in  the  metallurgical  history  of  ytterbium  foils 

can  be  a major  source  of  error.  Some  agreement  on  standardization  would 

be  helpful.  In  no  case,  however,  should  any  foil  but  cold-worked  (as-rolled 

or  unannealed)  be  used  as  gage  element  material, 

* The  melting  point  of  ytterbium  is  824°C. 

+ Obtained  from  American  Rare  Earth  and  Foil,  Inc.  (AREF) . 
tt  This  agreement  is  further  corroborated  by  Spataro's  (1972)  data  on 
0 . 0025-cm  (1-mil)  AREF  foil  up  to  4.7  kbar. 
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1-7.5  Double-Shock  Experiments— Shots  S3,  S5,  S7 , S9,  S29,  and  S30 


were  done  to  produce  peak  stress  loading  data  and  some  information  on 

whether  or  not  the  loading  calibration  depended  on  the  path  by  which  the 

gage  reached  peak  stress.  These  were  the  double-shock  experiments  described 

in  Section  4. 4. 4.  The  first  four  shots  we.s  designed  so  that  the  gage 

was  loaded,  released  to  zero  stress,  and  then  reloaded  to  a lower  stress. 

The  reproducibility  of  the  AH/H  values  associated  with  the  first  stress 

o 

peak  was  excellent  in  all  of  these  shots,  as  can  be  seen  from  the  data 
presented  in  Table  4.  The  poorest  accuracy  was  shown  by  gage  2 in  Shot 
S6,  where  the  deviation  from  the  expected  value  of  stress  was  7.9^,  but 
this  represented  an  error  of  only  0.01  kbar  at  this  stress  level. 

All  of  the  first-peak  data  from  this  series  were  accepted,  but, 
even  though  the  second-peak  data  looked  reasonable,  they  were  rejected 
tor  the  main  calibration  curves  for  the  following  reasons.  When  analyzing 
the  reload  data,  we  checked  back  on  previous  work  (Ginsberg,  1971)  and 

* 

Note:  The  corrected  data  from  our  previous  work  (Ginsberg,  1971) 

is  shown  in  the  following  tabulation. 


Shot  No. 

AH/H 

o 

L'ncorrected  g 
kbar 

Corrected  ^ 
kbar 

S3 

0 .0196 

0.78 

0.84 

0.0229 

0.78 

o 

00 

SI 

0.0160 

0.77 

0.82 

S5 

0.0208 

0.79 

0.87 

0.0218 

0.79 

0.87 

S6 

0.0225 

0.90 

1.03 

0.0218 

0.90 

1.03 
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found  that  an  error  was  made  in  assigning  stress  values  to  the  measured 

<\K/R  data.  This  arose  beeause  of  an  error  in  the  impedance  value  used 
o 

for  aluminum,  a number  needed  for  the  impedance  match  calculation 
neeessary  to  find  the  stress  level  of  the  second  peak.  The  analysis 
was  based  on  an  assumption  that  the  aluminum  was  elastic  in  the  stress 
range  involved  (below  3 kbar)  and  so  the  longitudinal  instead  of  the 
bulk  impedanee  should  have  been  used.  When  our  old  data  were  corrected, 
we  found  that  the  c libration  curve  in  the  region  below  1 kbar  was 
shifted  considerably  to  the  left.  However,  having  no  shock  data  in 
this  stress  range  for  comparison,  we  could  not  ascertain  whether 
the  value  of  AR/R  associated  with  the  reloading  stress  was  valid 
for  first  loadings.  Our  data  from  this  project  help  only  insofar 
as  all  of  the  old  second  shock  data  show  a lower  piezoresistance 
coefficient  than  any  of  the  first  shock  values  obtained  in  the  present 
project,  but  this  is  not  conclusive  evidence,  because  there  is 
probably  a small  difference  in  the  stat'3  of  work  hardening  of  the 
0 . 0025-cm  (1-mil)  foil  used  previously  and  the  0.005-cm  (2-mil)  foil 
used  in  this  project.  Within  the  present  series  of  shots  (S3,  S5,  S7, 
and  S9) , the  results  are  not  conclusive  either,  but  only  because  of  a 
laek  of  data.  There  seems  to  be  no  indication  that  second  shock 
AR/R^  values  are  either  greater  or  less  than  first  shock  values  for 
a given  stress  level.  To  be  safe,  however,  all  second  shock  data 
were  excluded  from  the  calibration  eurves. 

Shots  S29  and  S3C  are  also  described  in  Section  4.4.4.  In  each 
shot  we  attempted  to  compare  directly  the  response  of  a pair  of  gages 
that  reached  a given  stress  through  a single  loading,  and  a pair  that 
reached  the  same  peak  stress  in  two  steps.  The  results  as  shown  in 
Table  4 suggest  that  there  is  a small  effect,  based  on  the  observation 
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that  among  the  data  points  accepted,  the  £K/K  level  reached  through 

o 

a single  shock  is  slightly  greater  in  each  case  than  the  level  reached 
after  a step  load.  However,  the  effect  is  not  large,  and  we  cannot  at 
this  time  suggest  any  corrections. 

1 .7.(3  Phase  Transition  Experiment — The  records  and  a plot  of  the 
data  from  the  phase-transition  experiment  (S20)  are  shown  in  Figures  15 
and  Hi.  The  ytterbium  grids  each  had  1000-./  resistors  in  one  of  the 
signal  legs  and  the  manganin  gage  had  a 50-  resistor  in  one  of  the 
signal  legs.  This  is  the  major  reason  for  the  observation  that  the 
voltage  signal  from  the  ytterbium  gages  was  only  three  times  greater 
than  that  from  the  manganin  gage.  Figure  10  plots  the  output  from  the 
ytterbium  gages  and  the  output  from  the  manganin  gage  multiplied  by  a 
1 actor  of  at)  so  that  the  scales  would  be  comparable.  The  measured  peak 
stress  reached  during  this  experiment  was  11.5  kbar  assuming  linear 
behavior  of  the  manganin  gage.  The  stress  level  calculated  from  the 
measured  impedance  of  the  soda-lime  glass  used  for  the  projectile  head 
and  target  was  44.7±0.2  kbar  assuming  linear  elastic  behavior  of  the 
glass.  This  result  indicated  that  the  glass  is  no  longer  behaving  as 
a linear  elastic  solid  in  this  stress  range  or  that  the  coefficient  of 
manganin  is  not  strictly  linear,  or  both.  Insofar  as  the  exact  peak 
stress  is  not  of  critical  importance,  we  will  accept  the  41.5  kbar  value. 

The  most  interesting  and  important  feature  of  Figures  45  and  46  is 
the  difference  in  the  shape  of  the  manganin  gage  trace  and  the  two 
ytterbium  gage  traces.  The  manganin  gage  records  the  actual  shape  of 
the  input  waveform,  showing  a rise  to  about  41  kbar,  a flat  top  lasting 
nearly  2 ,-sec,  and  the  release  back  to  zero  stress.  The  ytterbium  gages 
both  rise  with  the  manganin  gape  to  an  average  indicated  stress  of  about 

* Two  points  were  rejected  because  excessive  dispersion  and  flow 
interruption  by  the  gage  packages  prevented  the  waves  from 
equilibrating  at  a constant  stress  level  before  the  arrival  of 
release  or  recompression  waves. 
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33  kbar  and  then  drop  off  steeply  to  half  the  maximum  value  of  resistanee 
change.  The  decrease  then  becomes  more  l'adua  1 at  this  time  and  continues 
smoothly  until  the  arrival  of  the  release  wave  between  2 and  2.5  p,sec 
when  there  is  another  slope  change.  This  new  slope  apparently  represents 
a combination  of  the  effect  of  the  continuing  phase  change  (until  some 
stress  is  reached  at  which  it  stops)  and  the  normal  decrease  in  resistance 
of  the  untransformed  material.  The  bump  occurring  at  about  3.5  asec  is 
a fracture  signal  arising  from  a wave  reflection  from  new  free  surfaces 
caused  by  spallation  of  the  target  behind  the  gage  plane. 

The  difference  between  the  manganin  gage  trace  and  the  ytterbium 

gage  traces  clearly  shows  the  occurrence  of  the  FCC  — BCC  transition 

observed  under  hydrostatic  loading  at  approximately  39  kbar.  Hough 

extrapolation  of  the  portion  of  the  ytterbium  traces  before  the  arrival 

* 

ol  the  release  wave  to  £R/R  = 0 suggests  that  the  transition  would  be 

o 

complete  in  loss  than  10  ^sec , considerably  faster  than  one  would  expect 
from  the  experience  of  the  investigators  who  used  hydrostatic  pressure 
and  found  that  the  transformation  required  hours  to  complete.  It  would 
seem  that  either  the  nucleation  rate  or  the  growth  rate  of  the  BCC  phase 
is  dependent  on  the  state  of  strain  or  on  the  generation  and  concentration 
of  lattice  defects. 

At  one  time  it  was  believed  that  the  ICC  - BCC  transformation  involved 
a change  in  the  energy  band  structure  of  ytterbium  with  an  attendant 
shrinkage  of  each  atom,  allowing  the  transition  from  a close-packed  (FCC) 
to  a more  open  (BCC)  structure.  It  was  subsequently  shown,  however 
(McWhan  and  Jayaraman,  1963)  that  the  required  decrease  in  volume  could 
he  accomplished  in  this  manner  by  a change  in  the  effective  coordination 
number  of  the  BCC  phase  from  8 to  12  and  by  recalculating  the  atomic 


The  actual  resistance  ratio  of  the  BCC  phase  to  the  FCC  phase  is  0.8, 
but  this  exercise  is  done  simply  to  estimate  the  speed  of  the  reaction. 
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radii  based  on  the  new  coordination  number.  The  kinetics  of  the  transformation* 

also  surest  that  no  electronic  rearrangement  takes  nlace  because 

the  change  in  resistance  is  too  slow,  even  at  the  microsecond  level.  The 

main  (and  obvious)  conclusion  is  that  ytterbium  cannot  be  used  as  the 

active  element  in  any  dynamic  stress  transducer  above  the  stress  at 

which  t tie  phase  transformation  takes  place  because  the  metastable  ICC 

phase  cannot  be  regained  above  the  equilibrium  transition  point. 

1 • y • 7 Static  Compression  Experiments — Some  representative  results 
1 rom  our  static  compression  piezoresistance  experiments  were  shown  in 
figure  ;>().  figure  51  contains  some  measurements  done  by  Smith  et  nl.  (1972). 

It  shows  essentially  the  same  general  leatures-a  change  in  response  with 
stiess  cycling  and  a static  eoeilicient  in  reasonable  agreement  with  the 
measured  dynamic  coefficient.  In  figure  50,  the  system  was  held  at  stress 
for  1 min  at  each  point  except  for  once  during  each  loading  cycle,  when 
it  was  held  at  approximately  1 kbar  for  20  min.  No  change  in  resistance 
with  time  was  noted.  During  the  release  cycles,  we  also  held  once  for 

20  min  and  obtained  a change  from  £K  K 0.0313  to  Alt  R = 0.0251  or  a 

_ 1 ° o 

change  of  less  than  10  'I/msec,  far  less  than  is  of  interest  to  us.  The 

same  phenomenon  was  noted  on  return  to  zero  stress  3 hrs  later.  The 

rate  of  change  ol  resistance  was  of  the  same  order  as  that  calculated  above. 

* Note:  We  briefly  attempted  to  analyze  the  kinetics  of  the  transition  by, 

foi  example,  litting  the  AH  H (t)  data  to  the  Johnson-Mehl  equation 
(Burke,  1905)  which  is 

. „ n, 

y = 1 - exp  -kt 

where  \ is  the  Inaction  transformed,  k and  n are  constants  and  t is 
time.  AH  Hq  was  assumed  directly  proportional  to  the  fraction  trans- 
formed. We  did  not,  however,  find  good  agreement,  probably  because 
of  this  assumption  and  the  observed  very  fast  initial  trans-formation 
rate,  which  is  at  variance  with  the  Johnson-Mehl  analysis.  Since  it 
was  not  critical  to  the  main  part  of  the  program,  this  effort  was 


abandoned . 


Deformation  during  the  stress  cycling  results  in  the  creation  of  lattice 
defects.  We  believe  that  the  resistance  change  is  probably  a result  of 


thermally  induced  reduction  of  the  number  of  these  defects.  This  room 
temperature  annealing  is  a slow  process,  but  the  effects  are  definitely 
measurable  and  can  be  seen  in  other  metals  with  melting  points  around 
that  of  ytterbium,  even  in  cases  where  deformation-induced  defects  have 
a much  smaller  electrical  effect  than  in  the  case  of  ytterbium. 

Tensile  Strain  Experiments — The  results  obtained  in  the 
tensile  strain  experiments  help  considerably  in  the  attempt  to  under- 
stand the  pi ozoros i s t anco  ol  ytterbium  and,  in  a qualitative  way,  some 
of  the  more  dif ficult— to— interpret  features  of  some  stress  gage  records, 
however,  the  results  of  these  experiments  cannot  be  uncritically  applied 
to  the  interpretation  of  stress  gage  records  because  some  important 
questions  are  still  unanswered.  The  three  major  ones  are:  (1)  How  much 
tensile  strain  is  actually  coupled  into  the  field  gage,  and  is  it  truly 
biaxial?  (2)  Mow  much  bending  is  present?  (3)  What  are  the  effects  of, 
for  example,  Poisson  contraction  of  the  gage  package? 

The  shape  of  the  curves  (Figures  47,  48,  and  49)  obtained  in  these 

experiments  is  very  interesting.  If  one  assumes  a constant  resistivity 

and  uses  Equation  (8)  to  predict  resistance  changes  during  tensile  loading, 

the  calculation  will  show  that  the  resistance  change  is  always  positive 

because  the  length  of  the  specimen  increases  and  the  cross-sectional  area 

decreases.  If  one  then  goes  on  to  assume  perfect  plasticity  and  constancy 

of  volume  during  deformation  (with  no  resistivity  change),  the  prediction 

will  be  both  positive  aid  constant  (Keough,  1970),  The  initially  negative 

* 

slope  shown  in  the  present  data  on  ytterbium  can  be  accounted  for  only  if  the 
resistivity  change  is  negative  and  large  enough  to  overshadow  the  geometry 
effects.  This  is  reasonable  when  we  consider  that  all  the  initial 
deformation  is  elastic  and  so  the  major  effect  of  the  deformation  is  a 

* 

Nickel  also  has  a negative  elastoresistance  coefficient. 
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change  in  interatomic  distance*  We  know  that  the  major  change  in  resistance 
during  compression  is  due  to  resistivity  changes  (the  magnitude  of  the 
effect  cannot  be  accounted  for  in  any  other  manner)  and  that  in  one- 
dimensional strain  the  resistivity  change  reflects  a large  change  in 
interatomic  distance  and  an  accompanying  change  in  band  structure.  It 
is  not  surprising  that  a tensile  strain  would  reverse  this  effect  and 
th it  deformation  that  affects  the  band  structure  in  the  opposite  way 
from  compression  should  cause  negative  resistivity  changes.  Our  measured 
tensile  yield  strain  of  0.005-cm  ytterbium  foil  is  approximately  0.27%. 

At  this  strain  we  would  expect  to  see  a change  in  the  slope  of  the 
F-  R/R  plot  because  the  deformation  is  no  longer  purely  elastic  and  slip 
by  dislocation  movement  is  taking  place.  If  the  material  work-hardens , 
we  should  expect  some  additional  elastic  deformation  with  accompanying 
negative  resistivity  change,  and  competition  between  the  two  processes. 
Inspection  of  Figure  17  suggests  that  a slope  change  at  about  0.27%  strain 
is  not  unreasonable,  given  the  uncertainty  in  our  data.  The  minimum  value  of 

at  somewhat  less  than  1.50%  strain  is  reached  when  the  positive  resistance 
change  arising  from  geometry  effects  and  from  lattice  defect  generation 
balances  out  any  further  negative  resistivity  changes,  and  the  subsequent 
positive  slope  continued  to  reflect  the  predominance  of  detect  and  geometry 
induced  changes.  Unloading  the  specimen  at  any  point  on  the  curve  caused 
the  resistance  to  drop  back  toward  zero  on  a line  of  slope  equal  to  the 
initial  slope — probably  a consequence  of  recovery  of  elastic  strain. 

Without  more  information  we  can  only  say  that  (1)  stretch  and  bending 
phenomena  are  probably  the  cause  of  any  small  negative  voltage  excursions 
on  field  gage  records  and  (2)  no  matter  what  the  answers  to  the  question 
r aised  above,  stretch  and  bending  effects  are  probably  negligible  at 
stresses  at  or  above  about  3 kbar.  A detailed  discussion  of  the  role  of 
our  tensile  data  in  checking  for  self-consistency  between  shock  and  hydro- 
•stratic  data  and  determination  of  the  state  of  gage  strain  is  presented 
in  Appendix  A. 
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'•7-9  General  Comments --Our  characterization  experiments  consisted 
°f  r°SiStlvity*  temperature  coefficient,  and  static  mechanical  property 
measurements,  as  well  as  X-ray  crystal  structure  determination  and  an 
attempt  at  hardness  measurements. 

The  measured  resistivity  of  ytterbium  was  found  to  differ  by  as 
much  as  a factor  of  three  from  some  previously  reported  values  for  heavily 
cold-worked  material  (Williams,  1970)  and  from  the  annealed  material  studied 
as  part  of  this  program.  In  general,  materials  with  lower  ambient 
pressure  resistivity  tended  to  have  a greater  sensitivity  to  stress  than 
materials  of  higher  ambient  pressure  resistivity.  This  suggests  that  the 
response  of  ytterbium  is  very  dependent  on  lattice  defect  concentration, 
and  especially  on  the  number  of  defects  actually  produced  during  the 
dynamic  loading  process. 

When  the  temperature  coefficient  of  resistance  of  our  foils  was 
measured,  we  found  that  the  annealed  AREF  foil  had  a temperature  coefficient  of 
2. 95  x 10  V.  - C , the  value  of  the  KC  foil  was  1.35  x lo“3y:j-°Ct  and 
the  as-rolled  AREF  foil  had  the  lowest  coefficient,  0.51  x 10"3ft/>OC. 

Our  experience  with  the  difference  in  chemical  composition  of  random 
samples  of  KC  foil  taken  from  the  same  lot  suggests  that  an  assumption 
of  compositional  invariance  for  the  AREF  foil  might  not  be  wise,  so  the 
difference  in  coefficients  among  the  KC  and  AREF  foil  and  between  the  two 
samples  of  AUEf  foil  might  be  due  to  either  chemistry  or  metallurgical 
history  effects.  In  our  previous  work  a coefficient  of  about  0.6  x lo"3 

V>°C  WUS  m0aSllred  lor  saraPl<-’s  ^ 0.0025  cm  (1-mil)  foil,  comparable 

to  the  AREF  as-rolled  0.005  cm  (2-mil)  foil,  so  it  seems  obvious  that  this 
property  cannot  be  predicted  from  chemical  analysis  or  knowledge  of  the 
metallurgical  history  of  the  foil.  The  importance  of  this  measurement 
derives  from  the  possibility  of  errors  in  stress  measurement  arising 
from  temperature  effects  in  the  gage  materia 1-the  most  desirable  gage 
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material  would  show  the  least  effect  of  temperature  on  piezoresistance 
coetficient,  and  would  probably  have  a low  temperature  coefficient  01 
resistance,  although  this  is  not  strictly  necessary. 

Li 1 ley  and  Stephens  (1971)  studied  the  effect  of  temperature  on 
piezoresistance  under  hydrostatic  pressure  and  found  that  raising  the 
temperature  decreased  the  sensitivity  under  those  conditions.  In  our 
experiments  the  stress  pulse  reached  t’’e  gage  after  current  had  pulsed 
through  it  for  about  10-15  asec.  The  field  gages  are  pulsed  for 
considerably  longer  times,  but  often  with  lower  currents.  Our  experience 
has  been  that  the  maximum  measured  temperature  rise  in  a 6-:.  encapsulated 
field  gage  made  from  0.005-cm  (2-mil)  foil  is  less  than  20°C,  using  a 
coeiiicient  ol  0.0  x 10  ~J .r  C.  This  rise  occurred  at  10  msec  with  a 
voltage  of  300  volts  across  the  bridge  of  which  the  gage  was  a part. 

At  200  and  100  volts,  the  temperature  rise  was  about  10°C  and  2°C , respectively. 
These  changes  in  temperature  should  have  a very  small  effect  on  the 
response  to  stress,  and  we  also  believe  that  there  should  be  no  temperature 
problem  in  transferring  our  laboratory  calibration  data  to  the  interpretation 
of  field  records. 

The  mechanical  property  measurements  were  done  because  we  have 

attributed  some  of  the  behavior  seen  in  ytterbium  to  yield  point  effects, 

especially  the  knee  in  the  ^ versus  'll  H curve  between  1 and  1.5  kbar 

o ' 

and  the  lack  of  zero  stress  residual  resistance  at  these  low  stresses. 

Moreover,  we  expected  that  the  tensile  results  would  facilitate  under- 
standing of  the  electrical  behavior  under  tension.  The  results  of  static 
compression  tests  (Table  3)  showed  ytterbium  ingot  to  have  a static 
compressive  yield  strength  of  about  0.3  kbar.  Cold-worked  foil  would 
have  a higher  static  yield  stress,  and  1 to  1.5  kbar  s>fcms  to  be  a 
reasonable  estimate  of  the  dynamic  compressive  yield  of  our  as-rolled 
material . 
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fj.  KKCOMMKNIJAIIONS  FOR  »•  L'Tl'UK  WORK 

\1 though  a considerable  amount  of  work  was  done  during  this  program, 
three  major  areas  remain  in  need  of  further  investigati  n.  These  are: 

• Quality  control  improvement 

• I’nders  t and  i ng  the  effects  of  complex  strain  states, 
such  as  those  encountered  in  the  field  and  reloading, 
on  gage  response 

• Data  treatment  and  development  of  a standard  method  of 
converting  gage  records  to  stress-time  histories 

Improvement  in  finality  control  methods  for  gages  and  gage  material 
would  lead  to  increased  reproducibility  and  reliability.  Some  relatively 
simple  test  that  could  he  used  to  identify  potentially  unsuitable  material 
would  he  the  goal  ol  such  a program. 

A better  understanding  ol  the  effects  of  complex  strain  states  such 
as  those  present  in  divergent  flow  fields  and  reloading  would  help  in 
the  analysis  and  interpretation  of  stress-time  histories  obtained  with 
ytterbium  gages. 

Hie  data  generated  during  this  program  is  extensive.  It  should  be 
examined  and  analyzed  in  more  detail,  and  a best  way  of  using  the  data 
lor  the  interpretation  of  gage  records  should  be  designed.  This  should 
then  become  the  standard  method  of  converting  gage  records  to  stress-time 
histories . 


116 


REFERENCES 


r 


L.  M.  Barker  and  R.  K.  doll  on bach,  1970,  Journal  of  Applied  Physics, 

V ol , ‘11,  p.  10, 

L.  Barsis,  E.  WilMams,  and  (.  Skoog,  1971,  Journal  of  Applied  Physics, 

Vol.  *11,  p.  5155. 

S.  11.  Brown,  1970,  "Development  of  Piezoresi st ive  Stress  History  Gages," 
Lawrence  Radiation  Laboratory  Livermore,  California,  Memo  DLT  70- 
*18. 

E.  Bucher,  P.  11.  Schmidt,  A.  Jayaraman,  K.  Andres,  J.  P.  Maita,  K.  Nassau, 
and  P,  D.  Dernier,  1970,  Physical  Review  B,  Vol.  2,  No.  10,  pp.  3911- 
3917. 

J.  Burke,  1965,  The  Kinetics  of  Phase  Transformations  in  Metals, 

Pergamon  Press,  Oxford,  p.  45. 

A.  It.  Daane,  1961,  in  The  Rare  Earths,  F.  H.  Spedding  and  A.  H.  Daane, 
od s . , John  Wiley  & "ons,  Inc.,  New  York,  p.  177. 

G.  E.  Dieter,  Jr.,  1961,  Mechanical  Metallurgy,  McGraw-Hill  Book  Company, 

New  York,  p.  480. 

D.  C.  Erlich,  D.  C.  Wooten,  and  R.  C.  Crewdson,  1971,  Journal  of  Applied 
Physics , Vol.  42,  No.  13,  p.  5495. 

M.  J.  Ginsberg,  1971,  "Calibration  and  Characterization  of  Ytterbium 

Stress  Transducers,"  Final  Report,  Contract  DASA01-69-C-0014, 

Stanford  Research  Institute,  Menlo  Park,  California, 

D.  E.  Grady,  et  al.,  1973,  "j_n  Situ  Constitutive  Relation  of  Rocks," 

Draft  Final  Report,  Contract  DNA001-72-C-0056,  Stanford  Research 
Institute,  Menlo  Park,  California  (January  1973). 

II.  T.  Hall,  J.  D.  Barnett,  and  L.  Merrill,  1963,  Science,  Vol.  139, 
p.  111. 

H.  T.  Hall  and  L.  Merrill,  1963,  Inorganic  Chemistry,  Vol.  2,  p.  618. 

J.  R.  Hearst,  1974,  Algorithms  for  Reducing  Ytterbium  Stress  Gage  Data, 

Lawrence  Livermore  Laboratory  UCID-16491. 

C.  M.  Hurd  and  J.  E.  A.  Alderson,  1973,  Solid  State  Communications,  Vol.  12, 
pp.  375-377. 


117 


KEFEKKNCKS  (Continued) 


* • X.  Kaysor,  1970,  Physical  Review  Letters,  Vol . 25,  No.  10,  pp.  662-664. 

F.  X.  Kaysur  and  S.  I).  Soderquist,  1969,  Scripta  Metal lurgica , Vol.  2, 
p.  259. 

0.  0.  Keough,  1967,  "Piezoresisti ve  Stress-Time-  Transducer  Development 
and  Granite  Adiabat  Measurements  for  Project  Pile  Driver,"  Final 
Keport,  Contract  DA-49 -146-X/.-28,  Stanford  Research  Institute, 

Menlo  / ik,  California  (February  1967). 

D.  D.  Keough,  1968,  Procedure  for  Fabrication  and  Operation  of  Manganin 
Shock  Pressure  Gages,"  Final  Report,  Contract  AF29(601)-68-C-0038, 
Technical  Report  No.  AFWL-TR-68-57 , Stanford  Research  Institute, 

Menlo  Park,  California  (August  1968). 

D.  D.  Keough,  1970,  "Development  of  a Iligh-Sensitivi ty  Piezoresisti ve 

Shock  Transducer  for  the  Low  Kilobar  Range,"  Final  Report,  Cont  'act 
DASA01-69-C-0014 , Stanford  Research  Institute,  Menlo  Park,  California 
(March  1970). 

D.  D.  Keough,  C.  W.  Smith,  and  M.  Cowperthwai te , 1971,  "Constitutive 
Relations  from  In  Situ  Lagrangian  Measurements  of  Stress  and 
Particle  Velocity,"  Interim  Report,  Phase  I,  Contract  DASA-01-70-C- 
0098,  Stanford  Research  Institute,  Menlo  Park,  California  (January  27, 
1971). 

D.  D.  Keough,  R.  F.  Williams,  and  R.  A.  Pasternak,  1968,  "Piezoresisti ve 
Shock  Transducer  Study  Using  Lithium,"  Final  Report,  Contract  18- 
5816  for  Sandia  Corporation,  Stanford  Research  Institute,  Menlo  Park, 
California  (August  12,  1968). 

D.  D.  Keough  and  J.  Y.  Wong,  1970,  Journal  of  Applied  Physics,  Vol.  41, 
p.  8. 

R.  W.  Keyes,  1960,  "The  Effects  of  Elastic  Deformation  on  the  Electrical 
Conductivity  of  Semiconductors,"  in  Solid  State  Physics,  Vol.  11, 
Academic  Press,  New  York,  pp.  149-221. 

118 


V 


REFERENCES  (Continued) 


h.  M.  Lilley  and  D.  R.  Stephens,  1971,  TID-4500  (l’C-4),  "Elect, -leal 

Resistance  of  \tterbium  as  a Function  of  Temperature  and  Pressure," 
Lawrence  Radiation  Laboratory,  Livermore,  California. 

1).  13.  MacDonald,  1973,  "Piezoresistance  Effects  in  Calcium  and  Ytterbium," 
Ph.l).  Thesis,  Syracuse  University,  New  York,  1971;  Uj  iversity 
Microfilms,  Ann  Arbor,  Michigan,  1973. 

D.  G.  McWhan  and  A.  Javaraman,  1963,  Applied  Physics  Lettirs.  Vol.  3, 
p.  129. 

C.  F.  Petersen  and  D.  C.  Erlich,  1972,  "Dynamic  Properties  of  Rock 

Required  for  Prediction  Calculation,  Draft  Final  Report,  Contract 

DASA01-71-C-0094 , Stanford  Research  Institute,  Menlo  Park,  California 
(October  1972). 

J.  T.  Rosenberg  and  M.  J.  Ginsberg,  1972,  Bulletin  of  the  American 
Physical  Society,  Series  II,  Vol.  17,  No.  11,  p.  1099. 

L.  Seaman,  al.,  1969,  "Classification  of  Materials  by  Shock  Properties," 
Final  Report,  Contract  I29601-67-C-0087 , Technical  Report  No. 
AFWL-Th-69-96,  Stanford  Research  Institute,  Menlo  Park,  California 
(November  1969). 

C.  S.  Smith,  1958,  Macroscopic  Symmetry  and  Properties  of  Crystals," 
in  Solid  State  Physics,  Vol.  6,  Academic  Press,  New  York,  p,  175. 

C.  W.  Smith,  1972,  "Stress  Wave  Measurements  on  the  Mine  Dust  Explosive 
Test,"  Final  Report,  Contract  DNA001-72-C-0235 , Stanford  Research 
Institute,  Menlo  Park,  California  (August  1972). 

C.  W.  Smith,  D.  E.  Grady,  L.  Seaman,  and  C.  F.  Petersen,  1972,  "Constitutive 
Relations  from  jtn  Situ  Lagrangian  Measurements  of  Stress  aid  Particle 
Velocity,"  Final  Report,  Contract  DASA01-70-C-0098,  Stanford 
Research  Institute,  Menlo  Park,  California  (January  1972), 


119 


REFERENCES  (Concl’tded) 


P.  C.  Souers  and  G.  Jura,  1963,  Science,  Vol . 140,  p.  481. 

S.  J.  Spataro,  1972,  Shock  Studies  of  Ytterbium  Stress  Transducers," 
UOPKK  72-2/,  Lawrence  Livermore  Laboratory,  Livermore,  California. 

R.  A.  Stager  and  H.  G.  Drickamer,  1963,  Science,  Vol . 139,  p.  1281. 

D.  R.  Stephens,  1964,  Journal  of  the  Physics  and  Chemistry  of  Solids, 
Vol.  25,  p.  423. 

R.  D.  Stromberg  and  D.  R.  Stephens,  1964,  Journal  of  the  Physics  and 
Chemistry  of  Solids,  Vol.  25,  p.  1015. 

A.  Tanuma,  W.  R.  Datars,  11.  Doi , and  A.  Dunsworth,  1970,  Solid  State 
Communications,  Vol.  8,  p.  1107. 

R.  F.  Williams,  1970,  "Examination  of  Possible  Modifications  to  the 

Ytterbium  Shock  Wave  Transducer,"  N66001-70-C-0572 , Department  of 
the  Navy,  Naval  Undersea  RSd)  Center  (May  15,  1972). 


120 


Appendix  A 


THE  PIEZORESISTANCE  TENSOR  COEFFICIENTS  OF  YTTERBIUM 


Introduction 

Previous  attempts  to  correlate  stress-resistance  data  obtained 
with  piezoresi st ant  transducers  in  shock  wave  experiments  with 
hydrostatic  data  have  not  been  entirely  successful.  With  the  exception 
of  a paper  by  Barsis  et  al.  (1971),  explanations  have  centered  on  the 
differing  changes  in  geometry  suffered  under  uniaxial  shock  wave  con- 
ditions and  hydrostatic  conditions  and  differences  in  defect  production. 

An  often  ignored  factor  in  these  explanations  is  the  tensor  nature  of 
piezoresistance.  It  is  usually  assumed  that  the  piezoresistance  tensor 
is  spherical  and  hence  representable  by  a scalar  quantity.  Barsis  et  al, 
(1971)  considered  the  tensor  property  in  a paper  on  manganin,  and 
MacDonald  (1973)  has  treated  the  tensor  aspects  of  piezoresistance  in 
some  detail. 

This  tensor  nature  is  important  in  ytterbium  and  can  assist  in 

explaining  the  different  behavior  of  dynamic  and  hydrostatic  data.  For 

an  isotropic  elastic  material,  the  fourth-rank  piezoresistive  tensor 

contains  two  independent  terms  TT  and  V . This  is  analogous  to  the 

XX  X 2 

elastic  stiffness  tensor  for  an  isotropic  material.  The  terms  TT  and 
TT  can  be  determined  from  two  independent  static  experiments,  and  is 

X Ct 

done  in  the  following  analysis  with  our  one-dimensional  tension  data  and 
with  hydrostatic  data  of  Lilley  and  Stephens  (1971).  With  knowledge  of 
ff  and  TT  the  one-dimensional  strain  coefficient  realized  in  shock 

XX  X ^ 

wave  experiments  can  be  calculated  in  the  elastic  region  of  the  ytterbium. 
In  the  plastic  region  of  the  material,  the  resistance  change  can  be 
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separated  into  a stress-  Induced  reversible  part  and  a damage-induced 

irreversible  part.  The  coefficient  for  the  reversible  part  can  be 

calculated  from  and  ff  by  assuming  that  material  response  above 

the  elastic  limit  is  perfectly  plastic.  We  ha'-e  compared  the  predicted 

stress-resistance  curve  for  a reasonable  HEL  in  ytterbium  with  the  stress- 

induced  part  of  our  shock  wave  data  and  with  Lilley  and  Stephen's  hydrostatic 

data  in  Figure  A-l.  The  tensor  coefficients  TT  and  it  are  con- 

11  12 

sistent  with  and  explain  the  one-dimensional  strain  coefficient  of 
approximately  0.04  (VQ/kbar  and  the  hydrostatic  coefficient  of  approxi- 
mately 0.06  Q'Q  kbar  (Lilley  and  Stephens,  1971).  The  following  analysis 
presents  the  tensor  theory  necessary  to  generate  the  model. 


The  piezoresistive  tensor  has  the  same  representation  and  satisfies 
llio  same  symmetry  conditions  as  the  piezooptical  tensor  (Smith,  1958). 
For  cubic  symmetry  or  isotropic  materials  the  resistivity  can  be  written 


O = o 6 +01T  (7 

Mij  Mo  ij  Vijki  k l 


(A-l) 


where  p is  the  strain-free  resistivity,  6 is  the  Kronecker  delta, 
° ij 

n,  . the  piezoresistive  tensor,  and  iy  „ the  stress  tensor.  We  are 
ijk£  k£ 

interested  in  the  stress-induced  change  in  resistance  and  we  can  write 
Equation  (A-l)  in  the  form 


ffijk£  ^k l 


(A-2) 


It  is  convenient  to  go  to  the  simplified  matrix  or  Brugger  notation 
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STRESS  — kbars 


FIGURE  A-1  PREDICTED  PIEZORESISTANCE  OF  YTTERBIUM  AND 
HYDROSTATIC  DATA 
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*1 

— ■ »tJ  O <*.J  - 1.6)  (A-3) 

For  an  isotropic  material  ( ^uch  as  polycrystalline  ytterbium)  this  takes 
the  form 


where  the  piezoresist ive  coefficients  satisfy  the  isotropy  condition 


rr. 


1 1 


" IT. 


12 


= 7T 


44 


The  fourth,  fifth,  and  sixth  terms  on  the  left-har  side  of  Equation  (A-4) 
correspond  to  a current  in  one  direction  due  to  an  applied  electric 
field  in  another  because  of  the  shear  stresses  in  the  material.  These 
are  not  of  interest  in  the  present  analysis  and  we  will  consider  only 
the  upper  submatrix  in  the  piezoresist ive  coefficient. 
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An  clement  of  the  transducer  is  shown  in  Figure  A-2  exhibiting  the 
present  coordinate  system.  Not*  that  the  resistivity  in  the  z.  direction 
is  of  ultimate  interest. 


Since  we  wish  to  relate  t ho  resistance  and  not  the  resist ivit\  to 
the  stress,  the  change  in  element  geometry  (area  and  length)  must  be 
accounted  lor  Tlii  s can  be  done  simple  In  generating  a-  additional 
matrix  to  account  lor  t h<’  geomet  ry  change  The  change  in  resistance 
is  related  to  the  change  in  resistivity  and  geometry  by 


All.  P -c  AA 

1111 


(A-(5) 


whet  e i corresponds  to  the  x,  y,  or  / direction  The  geometry  changes 
can  be  related  to  the  strain  For  instance 


--1  AA 

A . .x 

= e and  = p * v 

x .x  A \ z 

o o 


The  strain  can  be  related  to  the  stress  through 


e 

i.) 


U 

E ^kk 


6 

i.) 


where  U is  Poisson’s  ratio  and  E is  Young’s  modulus.  The  iinal 
re.su It  is 
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(A-7) 
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(A-8) 


Elastoresi stance  Tensor 

The  elastoresistanco  tensor  relates  the  change  in  resistivity  to 
the  strain  tensor  rather  than  the  stress  tensor  and  thus  provides  con- 
siderable simplification  when  the  elastic  state  of  the  material  is  more 
easily  represented  by  the  strain  field  such  as  a uniaxial  strain  experi- 
ment. In  tensor  notation,  the  elastoresi stance  tensor  provides  the 
modulus  relating  the  resistivity  and  strain  field. 

*ij 

1 = M e 

P i j kX  kl 

o 

It  can  be  related  to  the  piezoresistance  tensor  by  using 

4>. 

— ^ = 77  ° 

P ijkl  kl 

o 

and  the  st ress-st rain  relation 

CT  , = C „ e 

kl  k£mn  inn 

to  obtain 
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There  tore 


¥>.  . 

1 i 

= n C e 

p i jk/,  k-imn  inn 


m = rr  c 

ijk-t  ijmn  mnk/ 


For  isotropic  material  the  following  relations  result 


M = (A  + 2p)  ff  ♦ 2A  IT  ( A— 9 ) 

Mlv  = Aff  + 2(A  + j*)n  ( A-10) 


2 pIT 


1 t 


(A-ll) 


The  coefficients  A and  p jre  the  Lame7 const  ants. 

A geometry  matrix  relating  the  resistance  change  to  the  strain  is 
obtained  with  the  product  of  G in  Equation  (A-8)  and  the  elastic 
stiffness  matrix.  The  result  is  the  geometry  matrix 


G 

ij 


-1 

1 

-1 


(A-12) 


The  desired  result  is  the  change  in  resistance  related  to  the  strain 


AR 

= M e + G e (A-13) 

R ij  j ij  j 

o 

where 
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(A-ld) 


Is  provided  by  Equations  (A-9)  and  (A-10). 


Wo  can  now  calculate  the  coefficients  of  piezoresistance  or  elasto- 
resi stance  tensor  from  appropriate  static  data  and  with  this  predict 
the  dynamic  coefficient  or  at  least  establish  consistency  between  present 
static  and  dynamic  data.  The  first  step  is  accomplished  with  our 
present  uniaxial  tension-resistance  data  and  the  hydrostatic  pressure- 
resistance  data  of  Lillcy  and  Stephens  (1971).  The  second  step  follows 
by  assuming  that  the  ytterbium  transducer  material  mechanically  deforms 
in  an  elast ic-perfoctly  plastic  manner. 


The  one-dimensional  tension  data  used  in  the  calculation  of  it 
, 11 
ant  ir^  are  shown  in  Figure  A-3.  Ytterbium  gages  were  bonded  to 

aluminum  dog-bone  specimens  which  were  strained  in  one-dimensional 

tension.  Notice  that  the  tension  coefficient  is  initially  negative. 

The  elastic  piezoresistance  coefficient  (the  broken  line  in  the  figure) 

can  be  estimated  by  drawing  a line  through  the  data.  For  one-dimensional 

tension  CT?  * 0 and  a*  = = 0.  From  Equation  (A-7)  this  gives 


£iR 


z 


R 

o 


+ 


1 + 20 
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) 


<7 

Z 


(A-15 ) 
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UNLOADING  PATHS 


FIGURE  A-3  FRACTIONAL  CHANGE  IN  RESISTANCE  VERSUS  TENSILE  STRAIN 


Thus,  from  the  tension  data 


*11  + 


1 + 2l> 


-0.012 


Given  the  following  data  on  ytterbium, 


Ixnigi tudinal  sound  velocity 

Shear  sound  velocity 

Bulk  sound  velocity 

Density 

Bulk  modulus 

Young's  modulus 

Shear  modulus 

First  Lamey  constant 

Poisson's  ratio 


C = 
s 

C = 
o 

Po  = 
B = 

E = 

p = 

X = 
u = 


1 . 88  mm  psec 
1.02  mm /pace 
1.46  mm/psec 
6.97  g/cm^ 
148  kbar 
187  kbar 
72.5  kbar 
100  kbar 
0.290 


we  can  obtain  the  first  piezoresistive  coefficient  of  \rb,  which  is 


*11  = ■°,02() 


Under  hydrostatic  stress 
the  expression 


a 

x 


a = a 
y z 


-P  and  Equation  (A-7)  gives 


AR 

z 1 

= "(77  + 277  )P 

R 11  1 2 3B 

o 


(A-161 


Lilley  and  Stephens  have  measured  the  initial  hydrostatic  coefficient  of 
resistance  of  ytterbium.  Their  result  gives 
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FIGURE  A- 4 COMPRESSIBILITY  OF  YTTERBIUM  AND  C-7  EPOXY 


1 


-0.062 


"ith  t,lc  known  bulk  modulus  and  the  previously  calculated  value  for  n 
the  second  piezoresistive  coefficient  is  obtained 

ff12  ' -°019  • 


To  obtain  the  dynamic  coefficient  it  is  assumed  that  the  elastic 
strain  in  the  ytterbium  is  one-dimensional . This  assumption  is  not 
evident,  since  the  element  is  usually  mounted  in  an  epoxy  resin  and 
since  lateral  spreading  of  the  ytterbium  might  be  expected  in  the  softer 
epoxy.  However,  the  assumption  can  be  justified  by  noting  the  comparable 
compressibilities  of  ytterbium  and  epoxy.  tn  Figure  A-4  the  P-V  curve 
for  ytterbium  is  compared  with  that  for  C-7  epoxy.  Homalite  epoxy, 
which  we  used  when  obtaining  our  shock  data,  is  about  10%  denser  than 
C-7  epoxy.  This  would  further  improve  the  comparison. 

Assuming  one-dimensional  strain  during  the  elastic  shock  loading, 
Equation  (A-13)  gives 

aR 

~ = (M  - l)e 

R 12  x 

o 

or 

aR  M - 1 

z 12 

_ Q 

R X + 2p  x 

o 

Using  the  calculated  values  of  and  if  the  elastic  shock  co- 

efficient provides  the  relation 
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-iR 

~ = -0.039  a 

R x 

o 


If  elastic-perf eetly  plastic  behavior  is  assumed  above  the  elastic  limit, 
there  is  no  further  deformation  of  the  unit  cell  and  the  coefficient 
is  the  hydrostatic  coefficient  at  that  pressure  obtained  by  Lilley  and 
Stephens.  In  this  region  the  dynamic  curve  has  the  same  slope  as  the 


hydrostatic  curve  but  is  offset  at 
static  curve. 


a constant  value  above  the  hydro- 


The  calculated  coefficient  is  compared  with  the  experimental  data 
in  the  following  manner.  Shock-induced  resistance  change  has  been  shown 
to  consist  of  two  parts;  one  part  is  due  to  the  stress  state,  the  second 
is  due  to  the  introduction  of  damage  during  the  shock  process.  The  latter 
part  remains  as  residual  resistance  after  the  material  has  been  reduced 
to  zero  stress.  This  residual  resistance  was  measured  by  Ginsberg  (1971) 
and,  in  Figure  A-l,  the  residual  resistance  was  subtracted  from  the  total 
resistance  to  obtain  a plot  of  the  stress-induced  resistance  change. 

Also  shown  in  Figure  A-l  is  the  hydrostatic  curve  obtained  by  Lilley  and 
Stephens  on  an  ytterbium  gage  obtained  from  SRI.  The  predicted  shock 
curve  determined  from  static  data  is  also  shown.  The  value  used  for  the 
HEL  was  estimated  from  the  tension  data  in  Figure  A-3. 


The  agreement  is  reasonably  good.  This  simple  theory  using  the 
lowest  order  piczorcsistive  tensor  coupled  with  the  assumption  of 
elastic-plastic  material  response  appears  to  explain  the  fundamental 
difference  between  dynamic  shock  and  hydrostatic  piezoresistance  data. 
The  theory  also  servos  as  a caution  when  ytterbium  stress  transducers 
arc  used  in  different  configurations  in  the  low  stress  region.  The 
coefficient  can  differ  by  almost  a factor  of  two  depending  on  the  state 
of  strain  achieved  in  the  gage  element. 
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Appendix  B 


CALCULATING  £K/H  GIVEN  AV/V  FOli  FOLK-LEAD  STliESS  GAGES 

ThL“  constant  cl,rrent  circuit  configuration  for  obtaining  the 
resistance  changes  of  the  active  element  of  a four-lead  gage  is  shown 
in  Figure  B-l.  The  problem  is  to  relate  the  measured  voltage  changes 
across  11^  to  the  actual  resistance  changes  of  1^ , the  active  element 
<>f  the  gage.  The  solution  presented  here  is  a DC  analyses  which  accounts 
for  both  the  nonconstant  current  in  the  main  circuit  and  the  small 
finite  current  in  the  monitoring  circuit  and  is  based  on  the  treatment 
bv  Wilkinson  (in  Kcough,  1968). 

In  f.gu-e  B-l,  through  Hg  represent  those  gage  resistances  which 
varv  during  the  loading  cycle.  1^  is  the  resistance  of  the  active 
clement;  1^  and  l(.{  are  the  resistances  of  the  piezoresi stant  portions 
<>f  the  current  leads  actually  in  the  applied  stress  field.  (The  rest 
of  the  current-lead  resistance,  if  any,  is  constant  either  because  it  is 
nonpi ezoresi st ant  or  because  it  is  outside  the  area  being  loaded  and  it 

imped  into  n7  and  K^. ) Similarly  H ^ and  ltg  are  changing  portions 
of  the  signal-lead  resistances,  and  Kg  contains  the  constant  portions. 

V°  \o  ai*e  instant  resistances.  1^  is  the  ballast  resistance 
in  die  power  supply.  It?  and  1^  are  the  external  ballast  resistors 
Plus  the  Portions  of  the  current  lead  resistances  which  are  constant. 

'{9  18  tho  coftxlnl  sisnal  cable  termination  resistance  at  the  gage  plus 

the  constant  portions  of  the  signal  lead  resistances.  li  is  the 
resistance  across  which  the  voltage  is  being  monitored,  usually  either 
the  internal  resistance  of  the  oscilloscope  or  the  scope  resistance  in 
parallel  with  a 51-ohm  cable  termination  resistor  (in  which  case  the 
scope  resistance  can  be  neglected). 

for  this  circuit  the  voltage  V across  H is  given  by 
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V = 


R H E 
1 10 


(H  + R + It  f H + R + It  ) (It  + R f R + R + R ) - It 
1'J6781159  10  1 


(B-l) 


The  problem  is  to  calculate  Alt/ It  as  a function  of  AV/V  given  the 

° o 

initial  values  ol  It  to  R with  It  to  R varying  together  such  that 

± 1 5 


(Alt  /It  ) , (Ar  /it  ) 
1 i.0  .1  j,o 


i - 1»  ••••  5 and  j = 1,  . .,,5. 


The  subscript  o is  being  used  to  denote  initial  value  before  loading,  and 
the  subscript  on  the  relative  resistance  change  being  calculated  has  been 
dropped  since  all  the  relative  resistance  changes  are  assumed  equal  (but 
it  is  understood  that  the  resistance  of  interest  is  R^ , the  active  element) 
Note  that  since  only  relative  voltage  and  resistance  changes  are  being 
calculated,  h in  Equation  (B-l)  will  cancel  out  and  may  be  dropped. 

This  equation  may  be  solved  numerically  using  the  RESIST  program 

reproduced  in  BASIC  at  the  end  of  this  appendix.  The  input  data  to 

RESIST  consist  of  the  initial  values  of  the  10  resistances  plus  the 

largest  1^  to  be  calculated,  and  the  size  of  the  increments  by  which 

is  to  be  increased.  The  printout  gives  the  input  plus  AR/R  , Av/V 

o o’ 

and  (Alt/It  )/(Av/V  ). 

o o 

The  term  ( AR/IO /( Av/V^)  is  the  correction  factor  by  which  an  experi- 
mentally observed  AV/V  should  be  multiplied  to  obtain  the  true  AR/R  of 

° o 

the  gage  element.  It  is  an  increasing  function  of  the  observed  signal 

AV  Vo.  Note  that  lor  data  reduction  purposes,  a plot  of  the  correction 

factor  vs  Av/V^  is  a considerably  more  sensitive  method  for  obtaining 

R K than  the  more  direct  approach  of  plotting  AR/R  vs  AV/V  directly. 

o o 
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The  corrections  employed  in  the  calibration  data  presented  here 
ranged  from  about  1*  to  Hr.  When  very  large  changes  in  resistance  were 
anticipated,  the  value  of  If  was  increased  to  lower  the  correction 
factor.  This  also  lowers  the  signal,  but  it  was  more  than  sufficient 
for  our  purposes.  it  must  be  remembered  that  each  gage  user  should 
analyze  his  situation  with  respect  to  power  and  recording  and  see 
if  tli  i s correction  scheme  applies  to  his  system  before  using  it.  in 
particular,  this  analysis  does  not  apply  to  bridge-type  power  and 
recording  circuitry. 
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PROGRAM  RESIST 


20 
21 
22 
23 
2 A 

25 

26 
27 
26 

29 

30 

31 
50 
60 
70 
HO 
90 
92 
9 A 
96 

1 00 
1 10 
120 
1 30 
1 40 
1 50 
1 60 

1 70 
180 
190 
200 
20  5 

207 

208 

2 10 
220 
230 
500 
510 
999 

READY . 


REM  ••RESIST”  CALCULATES  THE  RATI  3 OF  THE  CHANGE  IN  VOLTAGE 
REM  TJ  THE  INITIAL  VOLTAGE  AS  A FUNCTION  OF  THE  CHANGF  IN 

REM  RESISTANCE  OF  R1  OVER  THE  INITIAL  VALUE  OF  R1  FOR 

REM  MANGANIN  RECORDING  CIRCUITS. 

REM  INPUT  - 

REM  A/B/C  - THE  INITIAL#  FINAL  AND  INCREMENT  OF  THE 

REM  RESISTANCE  R1 

KEM  RESISTANCES  R2#  R2#  HA,  R5#  H6#  R7»  R8#  R9#  AND  SM-K10) 

REM  OUTPUT  - 

REM  DELTA  R/RO  AND  DELTA  V/ VO  OF  Rl»  RATIO  OF 

REM  DELTA  R/RO/DELTA  V/VO 

KEM 

READ  A#  B#  C 
READ  R2#R3#R4#R5 
READ  R6»  R7#  K8#  R9#  S 1 
PRINT  ”R 1 -"/A 

PRINT  ”r<2#R3#R4#R5  - R2 # R 3#  R 4#  R 5 

PRINT  "R 6.*  7#  R8#  R9 # K 1 0 ="# R6# K7# R8# R9# SI 

PRINT 

PRINT  "DELTA  R/RO”»  "DEL  TA  V/VO"#"DF.LTA  R/RO/ DELTA  V/ VO" 

FOR  Ri-A  T*  A STEP  C 
LET  S2rR2cfM/A 
LET  S3*R3*R1/A 
LET  54=R4*R1/A 
LET  55=R5*R1/A 
LET  L = R 1 ♦ S 1 
M=R1»S2*S3»H6^R7»R8 
N = R 1 ♦ S 4*  S 5* R9*  S 1 
P~  R 1 ♦ R 1 
V-L/(<M*N)-P) 

IF  R 1 =A  THEN  500 
IF  V = VO  THEN  220 
*-  - < R 1 - A)  / A/(  < V-  VO  )/  VO) 

GJ  TJ  210 

PRINT  <R1~A)/A»CV~V0)/ VO • L 
NEXT  R 1 
G J TJ  999 
V0=  V 

GJ  TO  205 
END 
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